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STRATIFIED STEADY PERIODIC WATER WAVES * 

SAMUEL WALSH + 

Abstract. This paper considers two-dimensional stratified water waves propagating under the force of gravity 
over an impermeable flat bed and with a free surface. We prove the existence of a global continuum of classical 
solutions that are periodic and traveling. These waves, moreover, can exhibit large density variation, speed and 
- - - amplitude. 
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^ , 1. Introduction. Stratification is a common feature of ocean waves, where the presence 

of sahnity, in concert with external gravitational force, can produce substantial heterogeneity 
in the fluid. The pronounced effects that may accompany even a moderate density variation 
have earned stratified flows a great deal of scholarly attention, particularly in the realm of 
geophysical fluid dynamics. In this paper we develop a theory for stratified gravity waves 

P^ ' that are traveling and periodic. 

If we imagine a wave on the open ocean, past experience suggests that it may be regular 
in the following sense. First, it is essentially two-dimensional. That is, the motion will be 
identical along any line that runs parallel to a crest. Second, if we regard the wave in a 
coordinate system moving with some constant speed, it appears steady. Finally, the profile is 
periodic in the direction of motion, descending monotonicaUy from a single crest to a single 
trough once per period. 

#y^ ', We now formulate governing equations for waves of this form. Fix a Cartesian coordinate 

^ ' system so that the .x-axis points in the direction of propagation, and the y-axis is vertical. We 

'sj" I assume that the floor of the ocean is flat and occurs at y - -d. Let y - rjix, t) be the 

£^ ' free surface at the interface between the atmosphere and the fluid. We shall normalize rj by 

21 \ choosing the axes so that the free surface is oscillating around the line y = 0. As usual we let 

M = u(x, y, t) and v = v{x, y, t) denote the horizontal and vertical velocities, respectively. Let 

f^ \ p = p(x,y, f) > be the density. 

ir^ ' Incompressibility of the fluid is represented mathematically by the requirement that the 

f~^ I vector field (m, v) be divergence free for all time 

u^ + v^.-Q. (LI) 

rS [ Taking the fluid to be inviscid, conservation of mass implies that the density of a fluid particle 

j^ ■ remains constant as it follows the flow. This is expressed by the continuity equation 

p, + upx + vpy^O. (L2) 

Next, the conservation of momentum is described by Euler's equations 

{U, + UUr + VUy — i 

■ /, (1.3) 

where P - P(x, y, f) denotes the pressure and g is the gravitational constant. Here, of course, 
we assume that the only external force acting on the fluid is gravity. 



*This work was supported in part by NSF grant DMS-0405066 

^Division of Applied Mathematics, Brown University, 182 George Street, Providence, Rhode Island 02912 
(samuel.walshObrown . edu) 

1 



2 SAMUEL WALSH 

On the free surface, we must ensure that the pressure of the fluid matches the atmospheric 
pressure of the air above, that we shall denote Patm- Thus, 

P = Patm, ony^?](x,t). (1.4) 

The corresponding boundary condition for the vector field is motivated by the fact that fluid 
particles that reside on the free surface continue to do so as the flow develops. This observa- 
tion is manifested in the kinematic condition 

v-rjt + urjx, ony - rjixj). (1-5) 

Since we cannot have any fluid moving normal to the flat bed occurring at y = —d, we require 

v = 0, ony^-d. (1.6) 

Note that there is no accompanying condition on u because in the inviscid case we allow for 
slip, that is, nonzero horizontal velocity along solid boundaries. 

We seek traveling periodic wave solutions (m, v,p, P, rf) to (ll.lb -( ll.6b . More precisely, we 
take this to mean that, for fixed c > 0, the solution appears steady in time and periodic in 
the x-direction when observed in a frame that moves with constant speed c to the right. The 
vector field will thus take the form u - u{x - ct, y), v - v{x - ct, y), where each of these is 
L-periodic in the first coordinate. Likewise for the scalar quantities: p{x,y, t) - p(x - ct,y), 
P = P(x - ct,y), and 77 = //(x - ct), again with L-periodicity in the first coordinate. We 
therefore take moving coordinates 

(x — ct,y) I— > (x,y), 

which eliminates time dependency from the problem. In the moving frame (ll.ll )-( ll.3b become 

Uji + Vy = 

(u - c)px + vpv = 

(m - c)Ux + VMy = ^ •. ■ -* 

{u-c)v^ + VVy = -^-g 

throughout the fluid domain. Meanwhile, the reformulated boundary conditions are 

f V = (m - c)ri^ on y = rjix) 

I V = on 3; = -^ (1.8) 

i P = ^atm on y = riix) 

where u, v,p, P are taken to be functions of x and y, 77 is a function of x, and all of them are 
L-periodic in x. 

In the event that u - c somewhere in the fluid we say that stagnation has occurred, as 
in the moving frame the fluid appears to be stationary at that point. In order to avoid roll- 
up and other instability phenomena, we shall restrict our attention to the case where u < c 
throughout. 

Recall that we have chosen our axes so that 77 oscillates around the fine y - Q. Mathe- 
matically this equates to 



f 

Jo 



L 

T](x)dx^O. (1.9) 
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In effect, this couples the depth to the problem, so we need not treat li as a free parame- 
ter. The trade-off, as we shall see, is that this normalization will have significant technical 
consequences later. 

An often indispensable tool in the study of incompressible fluids is the stream function — 
especially in the case of two-dimensional flow in bounded domains. This is a function whose 
curl is the vector field (m, v) and thus, the gradient is orthogonal to the field at each point in 
the fluid. For a homogeneous fluid, a standard method is to then reformulate the problem for 
the stream function instead of the flow. What recommends this approach is the fact that the 
boundaries of the domain will be level sets of the stream function. Therefore, by adopting 
streamline coordinates, we can fix the domain and eliminate the free surface problem. 

However, for a heterogeneous fluid this generally proves insufficient, as it does not suf- 
ficiently capture the effects of stratification. Instead, we employ a more generalized object 
whose use was pioneered by Yih and Long, among others (cf. ||26l , IfTSl for example). Ob- 
serve that by conservation of mass and incompressibility, p is transported and the vector 
field is divergence free. Therefore we may introduce a (relative) pseudo-stream function 
ip - ip(x,y), defined uniquely up to a constant by: 

lA.v = - ^/pv, lAy = ^/p(u - c). 

Here we have the addition of a p term to the typical definition of the stream function for an 
incompressible fluid. This neatly captures the inertial effects of the heterogeneity of the flow 
(see the treatment in |26|, for example). The particular choice of y/p is merely to simplify 
algebraically what follows. 

It is a straightforward calculation to check that i// is indeed a (relative) stream function in 
the usual sense, i.e. its gradient is orthogonal to the vector field in the moving frame at each 
point in the fluid domain: 

(M — C)if/jc + Vt//y — 0. 

Moreover, (II.8I 1 implies that the free surface and flat bed are each level sets of i//. For defi- 
niteness we choose tfr = on the free boundary, so that tfr = po on y = -d, where po is the 
quantity 



= J yjp(x,y)[u(x,y)-c]dy. (1.10) 



PO 



To see that this is well-defined, that is, the integral on the right-hand side is independent of x, 
we calculate 



d.x ( ylp{x,y) [u{x, y) - c]) dy 

I f^"*' 

Vp(x, r]{x))v{x, T](x)) - I dy ( yfpv) dy = 0. 



We shall call po the (relatrve) pseudo-volumetric mass flux; it represents the amount of fluid 
flowing through a vertical line extending from the bed to the free surface and with respect to 
the transformed vector field ( -\fp{u-c), yfpv). The level sets of i/r will be called the streamlines 
of the flow. 

Since p is transported, it must be constant on the streamlines and hence, we may think of 
it as a function of xjj. Abusing notation we may let p : [po, 0] — > R^ be given such that 

p(x,y) ^ p{-ilj(x,y)) (1.11) 
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throughout the fluid. The reason for taking -i/' here will become apparent in the next sec- 
tion. When there is risk of confusion, we shall refer to the p occurring on the right-hand 
side above as the streamline density function. We shall focus our attention on the case where 
the density is nondecreasing as depth increases. This is entirely reasonable from a physical 
standpoint. Indeed, hydrodynamic stability requires that the depth be monotonically increas- 
ing with depth, making this assumption standard in the literature. The level set -if/ - po 
corresponds to the flat bed, and the set where -i/r - corresponds to the free surface. There- 
fore, we require that the streamline density function is nonincreasing as a function of -i//. 
By Bernoulli's theorem the quantity 



is constant along streamlines, that is 



E:^P + ^{(u-cf + v^)+gpy, 



(u - c)E^ + vEy = 0. 



This can be verified directly from ( 11.71 ). In the case of an inviscid fluid, E represents the 
energy of the fluid particle at {x,y). The first term on the right-hand side, P, gives the energy 
due to internal pressure. The second and third terms combined describe the kinetic energy, 
while the last is gravitational potential energy. When evaluated on the free surface, E is 
usually referred to as the hydraulic or total head of the fluid. 

Under the assumption that u < c throughout the fluid, and given the fact that E is constant 
along streamlines, there exists a function/? : [0, |/?ol] -^ R such that 

^{x,y)^m(x,y)). (1.12) 

dyj 

For want of a better name we shall refer to /? as the Bernoulli function corresponding to the 
flow. Physically it describes the variation of specific energy as a function of the streamlines. 
Define 



Jo 



B(p):^ /5i-s)ds 
Jo 

for po < p <0 and let B have minimum value Bmin- 

Let us briefly outline a few notational conventions. Let D,j denote the closure of the fluid 
domain 

D^ := {(x,y) eR^ :-d<y< 7](x)}. 

For any integer m > 1 and a e (0, 1), we say that a bounded domain D c R^ is C"^+" 
provided that each point in the boundary, denoted dD, is locally the graph of a C'"^" function. 
Furthermore, for fixed m> \ and a e (0, 1) we define the space C™^"(D) to consist of those 
functions / : D ^ R with Holder continuous derivatives (of exponent a) up to order m 
and that are L-periodic in the x-variable. Similarly, we shall take C'l.j{D) to be the space of 
m-times continuously differentiable functions which are L-periodic in x. 
Our main result is the following 

Theorem 1.1. Fix a wave speed c > 0, wavelength L > 0, and a relative mass flux po < 0. 
Fix any a e (0, 1), and let the functions /3 € C'^"([0, \po\]) and p € C'^"([/?o, 0]) be given such 
that the (L-B) condition holds (see Definition 13. 31 f or a precise statement). Also we assume 
the streamline density function p is nonincreasing. 
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Consider traveling solutions to the stratified water wave problem ( |l.l| )-( fL6] l of speed c, 
relative pseudo-mass flux po, Bernoulli function [5 and streamline density function p such that 
u < c throughout the fluid. There exists a connected set C of solutions (u, v,p, rj) in the space 
CI^"(D^) X C2+"(A^) X C2+«(A^) X Cl^^(R) with the following properties. 

(i) C contains a laminar flow (with aflat surface 77 = and all streamlines parallel to 
the bed) 

(ii) Along some sequence (u„, v„,p„, rjn) e C, either max-Q— m„ f c, min^j— «« i -co, or 
C contains more than one distinct laminar solution. 
Furthermore, each non-laminar flow (u, v,p, rj) &C is regular in the sense that: 

(i) M, V, p and rj each have period L in x; 

(ii) within each period the wave proflle rj has a single crest and trough; say the crest 
occurs at X — 0; 

(iii) u, p and rj are symmetric, v antisymmetric across the line x — 0; 

(iv) a water particle located at {x, y) with < x < | and y > -d has positive vertical 
velocity v > 0; and 

(v) 77' (x) < on (0, f ). 

A few remarks. The (L-B) condition above arises in the local theory as a means of guarantee- 
ing the existence of a minimizer to a particular Sturm-Liouville problem. In Lemma [331 we 
will show that a sufficient condition for (L-B) to hold is the following: 






gpiO)Po > j {7T^2B(/,) - 2B™„ + 2eo)'^' + 



(p - por((2B(p) - 2B™i„ + 2eo)'/^ + gp'(p))'jdp (L13) 

where 

eo:^mRx\{2g\\p'\\^ple^P«f\ {2g\\p'\\^f , (4||p'|U)', (8g\po\p(0)f']. (1.14) 

Though this is far from necessary, it has the advantage of being entirely explicit. 

In order to better situate this result in the larger context of geophysical fluid dynamics, 
let us give a quick comment on length-scales. As is evident from Theorem 11.11 our method 
allows us to treat mathematically waves with arbitrary speed and period. From a physical 
standpoint, however, we caution that our model ceases to be valid in certain scales. First, to 
ensure that the flow is incompressible, the Mach number must be far less than one. Also, 
in order to justify working within the inviscid regime, we must assume that the speed of the 
wave is substantially faster than the time scale for diffusion in the fluid (a situation typical 
of salinity in ocean wave, cf. |3| for a more careful derivation). Finally, to safely neglect 
Coriolis effects in ( 11.31 ), we need the rate of ambient rotation, O, to be of larger order than the 
ratio c/L. 

Finally, we note that stagnation will occur along some sequence in C unless serious 
pathologies occur. We shall prove in §7 that, if there is some {u„, v„, t]„} in C with 

liminfminMn J, -00, Hminfmin(c- u„) > 

then ?7„ -^ uniformly. That is, the fluid domain is pinching off' to nothing, while a current 
of arbitrarily fast leftward-moving particles develops. 

Now let us discuss the history of this problem. Physically, one of the most distinctive 
features of stratified flows is their ability to exhibit so-called internal waves. These are flows 
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in which the motion is essentially driven by a density gradient in the fluid rather than grav- 
itational force. Qualitatively, this can result in some highly counterintuitive behavior. For 
example, in the ocean one sometimes encounters "dead water", a phenomenon where a strat- 
ified wave is moving quite rapidly near the floor, but relatively slowly near the surface. A 
ship encountering such a wave would observe quiescent water, but experience a large amount 
of drag (cf. e.g. (E], for a lengthy discussion). Given such phenomena, previous mathemati- 
cal investigations of heterogeneous fluids have largely focused on the study of these internal 
waves, where the effects of density variation are most pronounced. The typical setup is to 
consider a two-dimension stratified fluid (or several layered homogeneous fluids) confined 
between two impermeable horizontal boundaries. Some substantial efforts in this genre in- 
clude ||2l[3][l5]|23l, among many others. 

The free surface problem has mainly been studied in the context of solitary waves, that 
is, waves that limit to a constant height at ±00. A scalar governing equation for the pseudo- 
stream function was developed by Long IfTSl in 1953, and later improved by Yih (cf. 1261 and 
the references therein) that greatly aided these efforts. The standard approach, for solitary 
waves, became to take the Long- Yih equation; assuming a given upstream density profile, 
one then work downstream to deduce properties of the wave. This is done, for example, in 
ll4l[T9l[T6landll24ll. 

On the other hand, traveling, periodic stratified waves have not received nearly as much 
treatment in the literature. The first substantial results are due to Dubreil-Jacotin in 1937 
ifTOl . generalizing her work on the homogeneous case in [9|. Remarkably, she had, even at 
that early date, already derived Long's equation (in fact. Long's equation is referred to as 
the Dubreil-Jacotin-Long equation in some circles). Dubreil-Jacotin was able to analyze 
the fixed boundary problem, i.e. existence of traveling periodic stratified waves between two 
horizontal plates. Linearizing around these solutions, she gave a system of integral equations 
governing small amplitude, traveling stratified gravity waves, from which she was able to 
produce a class of solutions of that type under certain assumptions. Yanowitch |25| later 
obtained similar results by using a variational argument with an alternate governing equation 
due to Love, again assuming small amplitude. In 1953, Ter-Krikorov |22|, by means of an 
asymptotic argument and Yih's equation, proved the existence of long, stationary stratified 
waves. This was in answer to a large body of results on the long wave problem that assumed 
the flow was homogeneous and potential. 

In this paper we take a new approach rooted in the work of Constantin and Strauss on 
the homogeneous case (cf. fS] |6l and others). In 2004, these authors were able to prove 
the existence of a large class of rotational (homogeneous) gravity waves that were regular 
in the sense we discussed earlier. This was done under very weak assumptions relating the 
volumetric mass flux to the strength of the vorticity f5l. By generalizing the methods of 
Constantin and Strauss to the heterogeneous case, we shall inherit many of the benefits of 
that paper. For instance, there will be no mathematical restrictions placed upon the wave 
speed or wavelength. Moreover, as the hypotheses of Theorem 11.11 make clear, we need only 
some inequality between the density variation, specific energy and volumetric mass flux in 
order to conclude existence of a global continuum of solutions. 

To emphasize the effects of heterogeneity, we hew closely to the organizational structure 
of Q. We begin, in §2, by using i/r to change variables and thereby fix the domain. We 
proceed to derive a scalar, nonlinear boundary value problem that describes the height above 
the flat bed in the new coordinates. As a consequence of (11.91 1. the reformulated problem 
takes the form of a quasi-linear elhptic differential operator added to an integro-differential 
operator. It is this second term that shafl require the utmost most care; we will show it 
completely accounts for the effects of the stratification. Given the long history of integral 
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equations in the study of heterogeneous waves, perhaps the presence of the integro-differential 
operator here can be viewed as natural. We emphasize, though, that one of the strengths of 
our method is that the governing equation is still "more-or-less" elliptic (in a sense that will 
become clear later) and thus, substantially easier to treat. 

In §3, we prove the existence of a 1 -parameter family of laminar flows. For these solu- 
tions the height equation reduces to an ordinary differential equation. However, because of 
the stratification term, this takes the form of a nonlinear boundary value problem with op- 
erator coefficients. Analyzing the linearized problem along the curve of laminar flows, we 
use a result of Crandall and Rabinowitz to prove that there exists a simple generalized eigen- 
value from which bifurcates a local curve of small amplitude solutions. As one might expect 
given the volume of work on small amplitude stratified waves, this proves substantiaUy more 
difficult than the corresponding results for the homogeneous case in fSl]. 

The analysis of §4 continues the local curve to a global bifurcation curve using a degree 
theoretic argument. We are thereby able to seamlessly treat both small and large amplitude 
solutions. However, this require strong a priori estimates in order to guarantee that the oper- 
ator has the necessary compactness for it to be admissible in the sense of degree theory. This 
will be achieved by first "freezing" the stratification term, then applying Schauder theory to 
the resulting uniformly elliptic operator. Finally, we estimate the frozen operator in order to 
return to the full problem. The main result of this section is a global bifurcation theorem in 
the form of an alternative, following in the footsteps of Rabinowitz 1201 . 

We devote §5 to investigating the nodal properties along the global curve. It is shown 
that C can only return to the curve of laminar flows in a certain interval of parameter values. 

In §6 we provide uniform bounds in the C^^"-norm along the global curve. These, with 
the results of §5, allows us to prove that Theorem 11.11 follows from the global bifurcation 
theorem. This is done in §7. 

2. Reformulation of the Problem. The main goal of this section is to reformulate prob- 
lem ( ll.7l i-( fTT8] ) so that the fluid domain D,, is transformed into a fixed domain, whose closure 
we shall denote R. As usual for two-dimensional incompressible flow, the main tool here 
will be the pseudo-stream function i// introduced in the previous section. Recall that we have 
i/r = on the free surface, ifr = -pQ on the flat bed, where po is the relative pseudo-mass flux 
and tfr is defined uniquely by the requirement that 

lA.v = - ^/pv, ipy = yjpiu - c). (2.1) 

In light of ( 12. Il l and (|l.ll )-( fT!2] ). the governing equations inside the fluid become 



-IpyilJxx + Ipx^xy = -Py - Pg 

whereas the boundary conditions ( 11.81 ) are 



in D^, (2.2) 



= -if/yT]x, on y = r](x) 

= ^atm, on y = j](x) (2.3) 

= 0, on y = -d. 

Recall we have that the quantity 

E^P + ^{(u-cf + v^)+gpy, (2.4) 

is constant along streamlines. In particular then, evaluating the above equation on the free 
surface i/^ = 0, we find 

IV.AI' + 2gp(x, 77(x)) (i](x) + d)^Q, on y = j](x) (2.5) 
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where the constant Q :- 2{E\^ - P^xm + gd)- Note that this Q gives roughly the energy density 
along the free surface of the fluid. Critically, however, we shall see that as a consequence of 
our normalization of rj, d is not a parameter for the problem. On the contrary, in all but the 
most trivial cases, d varies along C. In our analysis, therefore, we shall instead be viewing Q 
as parameterizing the continuum. 

We now wish to find from ( 11.7b and ( 12.41 ) a scalar PDE satisfied by i/f. In the case where 
the fluid is homogeneous, it is easy to show that this will take the form of a semihnear equa- 
tion: 

AiA = yiip) (2.6) 

for some function j : [0, \po\] -^ R. It is not hard to show that, in this scenario, y will describe 
the change in vorticity as a function of the streamlines. When one allows for density variation, 
however, one expects there must be an additional term accounting for the gravitational effects 
of the stratification. 

One can prove that the corresponding relation in the heterogeneous case takes the form 

dE dp 

di// di// 

This is known as Yih's equation or the Yih-Long equation. In this paper, it shall serve as 
our governing equation for iff. Recall that in the previous section we introduced the Bernoulli 
function /3 and streamline density function p. Rewriting the above expression we arrive at an 
equation of roughly the same form as ( 12.6b : 

P(if,)^Aif;-gyp'(-if,). (2.7) 

Were it not for the presence of y on the right-hand side, this would reduce to the homogeneous 
case — a fact that agrees with our intuitive notion that density stratification should reintroduce 
the depth into the problem as a serious consideration. We remark that comparing this to 
(12.6b we see that the final term on the right accounts precisely for the gravitational effects of 
stratification, the inertial effects having been captured in the choice of the ifr. 

With ( 12.7b in hand, we now make a change of variables to eliminate the free boundary. 
The new coordinates we will denote (q, p) where 

q^x, p^-i//(x,y). 

This scheme is sometimes referred to as semi-Lagrangian coordinates, in recognition of the 
fact that we are working, in some sense, halfway between streamline coordinates and the 
usual Lagrangian system (cf. II24I '). 

By means of a scaling argument, we may take L :- In. Then, under the transformation 

(x,y)i-> {q,p), 

the closed fluid domain D,, is mapped to the rectangle 

R := {{q,p) eR^ ■.0<q<27T, po<p < 0). 

The purpose of the minus sign is simply to flip the rectangle so that the free surface will 
correspond to the top of R, while the flat bed will mapped to the bottom. Given this, it will be 
convenient to put 

T := {(q,p) eR:p^Q}, B := {(q,p) e R : p ^ po). 
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Note that, in light of (11. lit and ( |1.12| i, we have that 

P^Pi-p), p^p(p). 

Moreover, the assumption that the streamUne density function is nonincreasing becomes 

Pp < 0. (2.8) 

Next, following the ideas of Dubreil-Jacotin (cf. ll9][T0l), define 

h{q,p):=y + d (2.9) 

which gives the height above the flat bottom on the streamline corresponding to p and at 
X = q. We calculate: 

1 K 

i/^, = - — , 'A.v = T^. (2.10) 

hp hp 

Note that this implies hp > 0, because we have stipulated that u < c throughout the fluid. The 
change of variables then gives: 

d, = d,-'^dp 

dy = hp dp (-2.11) 

dp = hpdy 

dq = d^ + hqdy. 

Using these expressions we can solve for u and v in (|2.10t to obtain 

1 h, 

u^c-^—, V--— ^. (2.12) 

^Jphp yfphp 

In order to formulate ( 12.71 1 in terms of h, we observe that 

Al/r = dlif/ + d]llf 

= -dx( y/pv) + dy{ yjpiu - c)) 
ih 






ip / \ lip / lip hp 

Hence, Yih's equation ( I2.7l l becomes the following 

- hlPi-p) = (1 + hl)hpp + hgqhl - Ihghphpq - g(h - d)hlpp (2.13) 

where we have used ( I2.9l l to write y - h- d. Recall, however, that we have normalized 77 so 
that it has mean zero. Taking the mean of ( 12. 9t along T, we obtain 

n-2n 

d^d(h)^ j- h(q,0)dq. (2.14) 

Jo 
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That is, the average depth d must be viewed as a linear operator acting on h. Namely, it 
is the average value of h over T. Where there is no risk of confusion, we shall suppress 
this dependency and simply write d. As we shall see, the addition of this integral term to the 
governing equations will be the single most significant departure from the homogeneous case, 
both technically and qualitatively. In recognition of this fact, we shall refer to the mapping 
h H^ -g{h - d{h))hl,pp as the stratification operator. Equation ( 12.131 1 will then take the form 
of a quasilinear elliptic operator added to the stratification operator 

Next consider the boundaries of the transformed domain. On the bed we must have by 
the definition of h that 

h = 0, onB. (2.15) 

If we were to reformulate the problem in terms of y not h, then the integral term would 
disappear in ( 12.131 1. but the boundary condition on B would become y = -d. We see then that, 
in the presence of density variation, the depth cannot simply be eliminated from the problem. 
Throughout the fluid we have by (12. It and (|2.12| i: 

Given this, the definition of Q given in (12.5b becomes the requirement 

l+hl + hl(2gph-Q)^0, onT. (2.16) 

Altogether, then, combining ( 12.13b . (12.15b and (12.16b we have that the fully reformulated 
problem is the following. Find (h, Q) e C^^"(R) x R satisfying 

( (1 + hl)hpp + hgghl - Ihqhphpci - g(h - d(h))hlpp = -h],l3{-p), po<p <0 

I l+hl + hl(2gph-Q)^0, p^O (2.17) 

[ h^O, p = po, 

and hp > 0. Here p e C'+°'([po,0];R+) and/? e C'+''([0, |pol];R) are given function with 
Pp<0. 

We now prove the equivalence of the height equation problem to the original Euler equa- 
tion formulation. Because the differential equations relating i// and h, (|2.2b -( |23T l, do not 
include p, this result follows more or less from the methods of the constant density case in 
||5j. Nonetheless, for clarity we recapitulate that argument here. 

Lemma 2.1. Problem (12.17b is equivalent to problem (Il.lb -( fir6] l 

Proof. The preceding development shows that ( ll.lb - (ll.6b implies ( 12.17b for some p, fi 
and Q. It remains only to prove the converse. Fix pa and let h be a solution to (12.17b of class 



(2.18) 



(2.19) 
throughout R. Also, we note that the free surface of the flow is given by rjix) - h(x, 0) - d{h). 



function y6, with hp > OinR. 








Define the C^i"(^) functions 








F(q,p): = 


1 


G(q,p) 


hq{q,p) 


hp(q,py 


hp(q,py 


Then, as hpg = h^p. 








Fg + FpG - GpF = 


-hpghf + hghpphf + 


hpghf - hpphqhf 
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Our first task is to recover the pseudo-stream function i}/. Fix xq e R and let ^ denote the 
solution of the ODE 

if/yixfi, y) = -F{xQ, -iA(xo, y)) (2.20) 

along with the initial condition t^(xo, rjixo)) - 0. By assumption, hp is bounded strictly away 
from zero on R. We may let (5 > be given such that F > 6. Thus t//(xo, ■) is increasing at 
a rate greater than 6 as y decreases. Thus ( 12.201 ) is solvable until ifr(xo,y) - -po, for some 
y. Then for each x e R, we may define i//{x, y) on some interval [y{x), //(x)] with y < rj. By 
uniqueness of solutions to ( 12.20b and the periodicity of h, it follows that ifr is periodic in x 
within its domain of definition. 

We claim that y{x) - —d{h) for all x e R. To prove this, let H{x,y) :- -G(x, -tfr(x,y)) 
for X e R, y e [^(x), t](x)]. Then by ( |ZT9] l and ( 1X201 1. 

Hy(x,y) = GpiAy = -GpF = -F, - FpG = -F, + FpH. 

On the other hand, by C' -dependence of ip on the parameter xq, we may differentiate to find 

i^x)y = {^y)x = -Fq + Fplff_,. 

Thus (/^v and H satisfy the same ODE. We have 4r(x, ri(x)) = by definition. Differentiating 
this relation yields 

(Ax = -i/'.v'?' = F(x, -il/(x,i](x)))?]'. 
Likewise, by definition of 77 we have 

77' (X) = hg, 

and thus 

H(x, Tj(x)) — —G(x, —ip(x, ri(x))) — F(x, — i/'(x, 77(x)))77'(x) 

by (12.18b . We have shown that H and i//^ have identical initial data at (x, ri(x)), hence by 
uniqueness we conclude 

4r,ix,y) = H(x,y) = -G(x, -ilr(x,y)), Vx e R, y e [y(x), //(x)]. (2.21) 

Now, C' -dependence of y on x enables us to differentiate the relation tfr(x,y(x)) - -po, 

dy 
dx 
By ( |Z20b and dZlTT i this implies 

G(x,-po)- F(x, -po) ^ = 0, Vx e R. 
dx 

But as we have seen, hg - 77', so that hg(q, po) = and therefore G(-, -po) = 0. Also, we have 

that F > (5 > 0, so we may conclude that y(x) = yo, a constant. 

Finally, we must show that yo = po- To do so we observe that h{Q,pQ) - 0, h(Q,Q) - 

77(0) + d{h). Then by ( 12.20b evaluated at xq = 0, we have 

m + d^h(0,Q)^flhpiO,p)dp = £^ 

^/(O) -dip(0,y) 



\l/x{x,y{x)) + tpy{x,y{x))-^ = 0, Vx e ! 






vo F{Q-4,{0,y)) 

- Coy 

= ^(O)-yo. 
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Thus yo - -d. 

It remains now to show that tfr constructed above constitutes a pseudo-stream function 
for a solution to (ll.ll l-( fTT6l l. We have already proved that i// - po on y - -d, and i/' = on 
y - ri{x). Also, 

f 2 + G^ = (/r2 + (/,2 ^ |V(/r|2, in ^ (2.22) 

by ( 12.18b and (12.21b . On the other hand, from the boundary condition ai p = in ( 12.16b . we 
find F^ + G^ - Q- 2gph on the free surface. But as q = x, and h - y + d,we may combine 
this with ( 12.22b to conclude that IVi/rp + 2gp{y + d) -0 ony - rjix). 

Finally, differentiating (12.20b and (12.21b and summing the two yields 



diji dtp 

' dy '' dx 



^^-Fp^+Gp^-G, 



- -FpF - GpG - Gq 

- hpphp^ + hpph^^h^^ - hpqhqhp^ + hqqh~^ - hqhpqKf. 

In light of (2.3), this becomes 

Aiff = -;6(iA) + gyp'i-ip), in R. (2.23) 

We now define u and v through ( 12.12b . and, by abuse of notation, take p{x, y) - p(-i/r(x, y)). 
Then, recalling the moving frame transformation and the definition of the pseudo-stream 
function, we have by (12.23b and the arguments of the preceding paragraph that, (m, v,p, rf) e 
C^+/'(D;) X C^e?(A^) X Clt"(D'n) X C'^,"(D^,) solves the original problem (Ol-dESll- □ 

3. Local Bifurcation. The goal of this section is to prove the existence of a local curve 
of small-amplitude solutions to ( ll.lb -( ll.6l ). The main product of our efforts will be the fol- 
lowing the theorem. 

Theorem 3.1. (Local Bifurcation) Let c > Q, pq < 0, a e (0, 1), and C^^" functions [i and 
p defined on [0, |/?ol]. [^O' 0] respectively be given. Consider the traveling solutions of speed 
c with relative mass flux pQ of the water wave problem (|l.lb -( fTT6l ) with Bernoulli function fi 
and streamline density function p such that u < c throughout the fluid. Ifji and p satisfy (L-B) 
(see Definition I3.3I ) and pp < 0, then there exists a C^ curve Che of small-amplitude solutions 
(u, v,p, 77) e CJ^^'iD^;) X CJ^^iD^) x C^,7(A^) x Cl^^iW^). The solution curve Cioc contains 
precisely one laminar flow (with rj = 0). 

3.1. Overview. Following the ideas set forth in the previous section, we shall work with 
the equivalent height equation formulation ( 12.17b . First we shall prove the existence of a 
1 -parameter family T of laminar solutions. We then show the existence of a local curve of 
non-laminar solutions bifurcating from a simple eigenvalue of the linearized problem along 
7~. The result will then follow from an application of the local bifurcation theory of Crandall 
and Rabinowitz. 

In several ways, it will be this section where the stratification term will be problematic. 
This will be immediately apparent in the next lemma, where its presence will make unattain- 
able an explicit formula for the laminar solutions. The consequences of this, at least from a 
technical standpoint, will cascade through the subsequent development, as we will be forced 
to make perturbation arguments along a curve of solutions we can only hope to represent 
implicitly. Though these concerns are purely mathematical in nature, they are not without 



STRATIFIED WATER WAVES 1 3 

physical analogue. Indeed, as we remarked earlier, one of the striking features of stratified 
flows is their propensity to exhibit significant internal waves — even in a small amplitude 
regime. 

3.2. Laminar Flow. Consider laminar flow solutions to the height equation ( I2.17l i. By 
laminar we mean parallel shear flows where the free surface is flat. Any such solution must 
then take the form H - H(p), with rj = 0. The height equation then reduces to an ODE with 
operator coefficients: 

Hpp - g(H - dmwlpp = -Hl/5(-p), po<P<0 (3.1) 

along with with the boundary conditions 

H^O, onp^po (3.2) 

and 

I + HJilgpH - Q) = 0, onp^O. (3.3) 

Owing mainly to the presence of the stratification term, ( |3.1t -( [T2] l constitutes a nonlinear, 
non-autonomous ODE boundary value problem which we cannot solve explicitly. Our ap- 
proach is motivated by the observation that if we can replace /? and p with functions of y, the 
resulting autonomous ODE becomes tractable. Indeed, any solution to the height equation 
should satisfy Hp > 0, hence the vertical variable y is a strictly increasing function of p, 
and so inverting the two is a valid change of variables. This is an adaptation of a technique 
commonly seen in the analysis of solitary waves in a channel, where it can be convenient to 
invert p and y on the profile at +oo (cf. 131 [161 l23ll and many others). As in the first section, 
define 



Jo 



B(p):^ /3{-s)ds, B,ni„:= min B{p)<0. 

P<^lPo,0] 

Lemma 3.2. (Laminar Flow) Suppose that the streamline density function p satisfies ( 12. 8t . 
Then there exists a 1 -parameter family of solutions H(-;X) to the laminar flow equation ( 13. It - 
( 1331 ) with Hp > 0, where < -2B,ni„ < A < Q. 

Proof. In accordance with ( |2.9l l, let F := H-d denote the vertical variable. By specifying 
the free surface 77 = 0, we force 7(0) = 0. Thus //(O) = d, and so by (13.31 1. 

1 + Ylilgpd - 2) = 0, on p = 0. (3.4) 

It therefore suffices to find Y satisfying 

Ypp + {l3{-p)-gYpi)Yl^Q, pQ<p<0 (3.5) 

along with the boundary conditions 

F = 0, on/7 = 0, (3.6) 

Y - -d, onp = po, (3.7) 

where d (viewed as a positive real number) and Y additionally satisfy (13.4b for some choice 
ofQ. 
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We now change variables. Put s := Y(p). Then from ( 13.51 ) we calculate 

. / dB dp\ 

(dB dp\ ds 

ds ds) dp 

Therefore we may rewrite (I3.5l l as 

(-\Yp^) +(f(Y)) =0, po<P<0, (3.8) 

where 



(dB dp 

F(Y) = {—-gY-f 

\ ds ds 



s=Y \ ds dp 



Y=s 



or equivalently 



ds \dp ds j ds 



For definiteness we take F(0) - 0, so that integrating we find 



r- 



F(Y)^B(p)+ I gAr)dr, (3.10) 

Iy ds 



and 



-0 

F(Y)^B(p)+ I gY(ry-^(r)dr. (3.11) 



Jp dp 



Note that by ( 12.81 1. the integrand on the right-hand side is nonnegative, hence F > B^in. 
Returning to ( 13.81 ) we integrate once to obtain 

Yp^ , ^ =. (3.12) 

^JA + 2F(Y) 

where A is a constant of integration with A > -IB^i^- Note that, by (13.101 ). F is actually an 
unknown. Using the fact that F{Q) - 0, ( 13.121 ) and ( 13.41 ) together imply 

Q-A 

-' ^ (3.13) 



2gp{Q) 
Inverting ( 13.12b and combining it with ( 13.9b we arrive at a first order system for p and F: 

i ^ = ^JA + 2F(s), 

(3.14) 



f = (^B'(p)-gsp'ip)^^jA + 2F(s), 



where by abuse of notation p and B here have been extended continuously so that their domain 
encompasses the entire real line. We have that both Y and p vanish when s = 0. Locally, then, 
we can solve this initial value problem. Let {F, p) - (F{s; A), p(s; A)) be such a solution for a 
fixed choice of A > -IB^in and let the maximum domain of definition be (s-{A), s+(A)). We 
must show that p attains the value po somewhere in (s_, 0). 
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Elementary theory of ordinary differential equations tells us that, if the domain of defi- 
nition is bounded from below, then the solution must blow-up as we approach the boundary. 
More precisely, we must have \F\^ + \p\^ — » oo as 5 — > i-, provided that s^ > -00. First 
consider the F equation in (I3.14l i. Rearranging terms and integrating in s yields 

F = ( J (B'ip) - grp'(p)) dr^ -A^ \F(s; A)\ < \A\ + {\mM + gs%'\\Jj . 

It follows that F is bounded for finite s. Next, recall that we have chosen A > -2Bniin , whence 

A + 2F > A + 2Bmi„ > 0. 
Puting e(A) :- A + 2Brnin , by ( 13.14b we have 

dp r- 

-f>Vi>0. (3.15) 

ds 

But this uniform lower bound on -j^ implies that, if i_ - -00, then p — > -00 as s — > .?_. On 
the other hand, if s- > -00, then the arguments of the previous paragraph guarantee that \F\ 
is bounded on (.^-jO], so again we conclude p — > -00 as i — > s-. In either event, we may 
choose d = d{A) > to be the unique number satisfying pi-d) - pQ. Inverting p (which is 
vaHd by ( 13.15b ) we get a function Y - Y(p; s) satisfying ( I3.5l )-( l3.7b . Finally, let Q - Q(A) be 
defined by the relation ( 13.13b . 

Q(A)^A + 2gp(0)d(A). (3.16) 

Thus for each A > -2B,^i„, H(p; A) := Y(p; A) + d(A) is a solution the laminar flow equations 
with Q = Q(A). D 



3.3. Eigenvalue Problem. Up to this point we have produced a family, T, of laminar 
solutions to the height equation parametrized by the variable A, which is drawn from a suitable 
range determined a priori by the given function/? and the value of po. We will now linearize 
the full height equation around H for fixed A by evaluating the Frechet derivative. Ultimately 
we show that the linearized problem has a simple eigenvalue at some A*. 

Fix e > and let h{q, p) - Hip) + emiq, p). Then from (12.17b 

(1 + e^n?^){Hpp + enipp) + eniqqiHp + enipf - 2e^mgmgp(Hp + enip) 

-g(H + em - d(H) - ed{m))(Hp + empfpp - -(Hp + empfl3{-p), 

for Po < p <Q. Differentiating in e and evaluating at e = yields the linearized equation: 
mpp + m,gHl - gpp{3mp(H - d{H))Hl + (m - d(m))Hf} = -3Hlmpf5(-p), (3.17) 

for Po < p <Q. 

Applying the same procedure to the two boundary conditions we find 

2gpmHl + 2{2gpH - Q)mpHp =0, on p = 

and 

m = 0, on p - pq. 
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We may simplify the nonlinear condition on T by noting that for all A, 



Also, from dTTOl i. 



^(0) = :p-7^ ^ 2gp{Q)H{Q) -Q^-A. (3.18) 



Hl-Yl^j^^^A-\ onp=0. (3.19) 



Combining these we find that the the linearized problem is the following: for fixed A > 
-2Binin, find m that is 27r-periodic in q satisfying 



' nipp + niq^Hl - gpp{3mp(H - d{H))Hl + (m- d(m))Hl) 

= -3Hlmpl3(-p), po<p<0 

m-0, P - Po 

gpm = A^^^nip, p -Q. 



(3.20) 



Observe that the only effect of the variable density on the boundary is in the addition of the 
constant p(0) on p - Q. This similarity with the constant density case will allow many of the 
same calculations to push through virtually unaltered. On the other hand, the stratification 
term in ( I3.17l i will present a significant technical barrier as it introduces both the nonlocal 
operator d, and a zero-th order term. 

From this section onward we shall consider only A in the range A > -2B„^in + eo, where 
€() is as in ( 11.141 ). As will be made clear later, the purpose of slightly decreasing our range 
of admissible A is essentially to ensure that Hp is bounded uniformly away from zero for A 
small. 

In order for bifurcation to occur, some control over the relative sizes of p, B, fi and po 
is necessary — even in the homogeneous case. With this in mind, we make the following 
definition. 

Definition 3.3. We say the pseudo-volumetric mass flux po, streamline density function p 
and Bernoulli function p collectively satisfy the Local Bifurcation Condition provided 



' -gp(m(Qf + C "p'^P' ^)Mp)'dp ' 



inf inf i j: '^^^- ■ \ < -1 (L-B) 

..-2z*„„„..„ ,er I j^j2^p^^Hp^^p;A) + 8Pp(p))dp J 

where ^ := {0 e H\{po,Q)) : * 0, (pipo) = 0) and H{-;A) is the solution to (l3T1l-(l33Tl 
given by Lemma [372] 

We shall restrict our attention to case where (L-B) is satisfied. The next lemma gives the 
precise motivation for this choice. 

Lemma 3.4. (Eigenvalue Problem) Assuming ( 12.81 1 and (L-B), there exists A* > -IBj^in + 
eo and a solution m(q, p) f of (I3.20l l that is even and 2n-periodic in q. 

Proof. We seek special solutions of the form m{q, p) - M(p) cos(kq) for some A: e Z^, as 
we require 27r-periodicity in q. Note that this implies 

p2;r 

d{m) = M(0) -T cos{kq)dq = 0. 
Jo 

Also, since solutions of this form will necessarily satisfy the ODE 

mgc, - -k^m. 
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we may rewrite the interior equation of (|3.20| i in self-adjoint form: 

{Hfmp)p + {{H-' + gppk-^)m^% = 0. 

In turn, this requires that M satisfy 

{H-''Mp)p = ieH-p' + gpp)M, po<p<0. (3.21) 

Any value of k will suffice, but for simplicity we focus on finding a solution where k = 1. 
To do so we approach ( |3.21| l as a Sturm-Liouville problem. Note that, as a consequence of 
( |1.14t . the term in parenthesis above is nonnegative for k = I. 
In that connection we consider the minimization problem: 

u = Li(A)^ in fRiS; A) (3.22) 

where the Rayleigh quotient K is 

-gp(0)4>(0f + /° Hpip; AT'Mpfdp 

n^; A) := — -. ^ (3.23) 

lj(p)^{Hp(p;A)-'+8Pp(p))dp 

and S^ :- {(f) e i/'((/?o,0)) : S 0, 4>(po) - 0). For each A, the function M that attains the 
minimum, p{A), will satisfy the boundary conditions of the linearized problem (13.20b as well 
as the ODE: 

{Hp'Mp}p = -p(A)(Hp' +gPp)M. 

The task is to find A* > -28^^ + ^o such that piA*) = -1. The corresponding M will be 
precisely the special solution we seek. 

By the continuity of p and in light of the (L-B), it suffices to prove that p{A) > -1 for A 
sufficiently large. By ( 12.81 ) we have thatpy, < 0, and we found in the previous section that this 
induces a uniform lower bound on F{-;A). That is, for any A > -2Bniin + £o, 

min F(Y(p; A); A) >B^n, 
p 

We may therefore choose A such that 

A > gpiO) + gWppWi + g''^ VmO)IKIU - 2B^i„ 



-'plloo 

and conclude 

Hp' ^(A + IFiYipW^ >{A + 2B,„i„)'/2 > ^||^^^||^ ^ ^/^^_ 

Then, again using the fact that pp < 0, we obtain 

Hp'+gPp> Vi^. 
Now letw e S^ be given and fix a /I > -2Bmin + £o- Then 

JrO pO 

(Hp'wl + (Hp' + gPpW)dp > Vip(0) (w^ + gp{0)w],)dp 
Pa Jpo 

> 2gp(0) I wwpdp 

•Jpo 

= gp(.0)w(Of. 
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Thus Hiw) > -1. As this holds for arbitrary w ^ y and all admissible A, we conclude 
li{A) > -1. D 

As we indicated in the first section, there is an explicit condition, (11.131 1, on p^, p and fi 
that implies (L-B). Moreover, for any choice of j8, p, taking po sufficiently small (and restrict- 
ing p and p to the decreased domain), will be enough to guarantee that this size condition will 
hold. Though the size condition is not necessary, it is still general enough to allow for a great 
variety of flows. 

Lemma 3.5. (Sufficiency of Size Condition) Ifpo, p and/3 satisfy the size condition (I1.13I I. 
then they satisfy (L-B). 

Proof. We must show that for some A > -IB^^^^ + eg, we have yu(/l) < -1 in the sense of 
( 13.221 1. In ( 13.151 ) we showed that -j- was bounded below uniformly in p. Therefore, 



(^ dp ^ 

-PQ - -p(-d) - -rds > y/eod. 
J-d as 



If eo = 0, then we are in the constant density case and this lemma has already been proved in 
Il5]. Otherwise, 

d<-^. (3.24) 

As Yp = Hp, we see by ( ITT2b 

,1/2 



H-p'{p- A)^(a + 2B(p) + 2 r gpp(r)YAr)dr^ 
<(A + 2B(p) + 2gd(A)\po\\\Pp\U''^ 



i^A + 2B(p) + eo) 



1/2 



where the last line comes from ( 13.241 ) and the definition of eo. Pairing this with the size 
condition (|1.13l l, we have that for A sufficiently near -2Bniin + £()'■ 



f {Hf +{p- pofiHp' + gpp))dp < f ((A + 2B(p) + eo)''^ 

J Pq J Pq ^ 

+(P - PofiiA + 2B(p) + eo)"^ + 8Pp)yp 

< gp(0)pl- 

Now take w = p - po & S^ and let A be such that the above inequality holds. Then, 

-gpiO)pl^CHp'dp 

^(w;i) = -g < -1. 

jp^ip - PoHH-^ + gpp)dp 

It follows that p{A) - infy R{-;A) < -1, as desired. Hence the Local Bifurcation Condition 
holds. D 

For convenience we denote G{p;A) :- 2F{Y{p;Ay,A). Thus, for instance, we have 
Hpip) - {A + G{p;A)Y^^^. The great difficulty that arises from the stratification term in 
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( 13.11 ), in large part stems from the fact that G depends on A. Given that, before proceeding 
further, we set ourselves to a technical task: We must determining precisely how the A+G{-;A) 
term varies along T. 

For the next several proofs we will denote differentiation with respect to /I by a dot. 

Lemma 3.6. (Sign of 1 + G) For A > — 2Z?,„,>, + eo, 1 + G is strictly positive. In fact, we 
have the following bound: -1/2 < G < 0. 

Proof. By rewriting ( 13.111 ), we have the following expression for G: 

G(p; A) = 2B(p) + 2 f gY(r, A)p„(r)dr. (3.25) 

Jp 

Since B is independent of A, and indeed the /t of G derives entirely from that of Y, we differ- 
entiate to find 

1 + Gip; ^) = 1 + 2 r gY(r; A)ppir)dr (3.26) 

Jp 

But note that from (3.6) we can calculate 

1 1 r° 

Yp = --{A + G(p)r'^Hl + G) ^ Yip) = 2 J ^^ + G(r)r"Hl + G(r))dr, (3.27) 
hence 

-0 rO 



l 1 *'■' 



1 + G(p) =1+1 ( 8P'(r)(A + G(s)r^'\l + G(s))dsdr. 



Now we reverse the order of integration to find: 

1 + G(p) = 1 + r r gpp{r)(A + G(s)r"\l + G(s))drds 

= 1+ f kip,s)(l+Gis))ds, (3.28) 

Jp 

where k(p, s) :- J' gpp(r)(A + G{s)y^^^dr. This is a linear Volterra integral equation of the 
second type, albeit a simple one. As A: is a continuous map on A := {(/?, s) : po < p < s <Qi] 
we have that 

1 + G(p) = 1+1 hp, s)ds, (3.29) 

Jp 

where 

klip, s) := k{p, s), k^ip, s) := ( k,„-i(p, r)k(r, s)dr, m>2 

Jp 

and the resolvent kernel k is given by 



kip,s) := 2_j^'"^P'^'>- 



m=\ 
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The fact that the series converges comes from the following estimate: 

\{s-p)"'-^ 



\k„,(p,s)\ <( sup \k(p,s)\y 

\ (p.s)eA I 



(m- 1)! 



From this it is clear that controlling the quantity in parentheses allows us to control the size 
k, and thereby the sign of 1 + G. However, as we are taking take A > -IB^i^i + eo, we have 

sup \kip, s)\ < 8\po\£[i^'^\\Pp\\oo, 

(p,s)eA 

since G > 2Bniin implies A + G > eo. Inserting this into the definition of ^ we find 

sup |^(;,,.)| < g\pQ\e,"%p\W-P < Ipoko'^'lKlUe'""!, 

(p,s)eA 

Integrating this expression and applying the definition of eo 

l+G(p)>l-ple^''^\\Pj,\U^P''^^^. 

Thus we have shown that 1 + G is strictly positive, as desired. 

Finally, to derive the upper bound we merely reexamine ( I3.27l i. That is, from ( 13.271 ) we 
see that F > 0, as we have shown 1 + G > 0. But then ( 13.26b implies G < 0. D 

With the sign of 1 + G established, we are now in a position to better understand the 
relation between Q and A set down by ( 13.16b . 

Corollary 3.7. (Convexity of Q) For all A > —2B^i„ + eo, we have that Q is a convex 
function of A with a unique minimum /lo. 

Proof. First, recall that by the definition, Q{A) -A- 2gp(Q)Y{pQ; A), where Y(-;A) is the 
corresponding solution to the laminar flow. Differentiating twice in A, we obtain 

- ^ ,r,, = Y(po;A)- 
2gp(0) 

It follows that to prove the corollary it suffices to show Y(po; A) <0 for all A > -2Bmm + eo- 
To do so we first note that by differentiating ( 13. 26b by A 

G(p;A)^2 f gY(r,A)p'(r)dr. (3.30) 

Jp 

On the other hand, differentiating ( 13.27b we find 

''^' 2jp (A + G(r,A))"^ ^Jp (i + G(rU))5/2 

Substituting (13.31b into (13.30b and exchanging the order of integration we find that for each 
fixed A, Y{-; A) solves the integral equation 



I JiP^i 
Jp 



Y{p)^i{p)+ j{p,r)Y{r)dr, pe[p^M (3.32) 

h 
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where 

J;; (/I + G(i))-^/^ 2 J,, (^ + G{f)yi^ 

Observe that f < and ||^||oo - -{(po)- Now this is a Voherra integral equation with solution 
given by 



I lip, I 



Yip) = {(p)+ j(p,r)mdr 
h 

where the resolvent kernel j is given in the same fashion as the previous lemma. From this 
representation formula we immediately derive the estimate 

FC/^oU) < (1 - IpoIIIJIUMpo). (3.33) 

As {{po) < 0, we are done if we can show that the quantity in parenthesis is positive. Unsur- 
prisingly, this is again a consequence of the definition of eo in (11.131 1. Again let A denote the 
triangular region A := {{p, r) : po < p < r < 0). Observe that 

.-3/2. 



sup \j{p,r)\ < fg gllPplloolpoL 

(p,r)eA 



and thus 



But by the definition of eo, this becomes 



sup |;>,r)|<|;,oko'''^IIPplUl"»l. 

(p,r)eA 



sup \]{p,r)\ < -— 

(p,r)£A AP0\ 

Inserting this last expression into (|3.33t proves convexity of Q (and the concavity of Y). 
To show the existence of a global minimum we note that 

Q(A)^l-2gp(0)npoU), 
thus Q > for A sufficiently large, since we have shown Y(po; A) < 0, for all A. D 

We now return to the linearized problem. With the previous lemma in hand we are pre- 
pared to prove the monotonicity of fi in some neighborhood of A* . This will in turn give the 
uniqueness of A* . For convenience, we denote a = a{p; A) := Hp(p; Ay\ 



Lemma 3.8. (Monotonicity of p) pis a strictly increasing function of A whenever iJ.{A) < 0. 

Proof. For each w e S^ and each A, define L^w := -{Hp^Wp}p. For each A, let w{p) - 
w(p; A) be the solution to the eigenvalue problem 

L,w = piHp' + gpp)w, wipo) = 0, Wp(0) = A-^'^gp(0)w(0) 

where p - p(A) is the minimum eigenvalue. As in the previous lemma, denote by a dot dififer- 
entiation with respect to A, and for notational convenience we drop the explicit dependence 
on A. Thus, a = ^. Likewise, 



■^Lw = -{^fl(l + G)wp + a'wp'^ . (3.34) 
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i( ■ ] d /I +G\ 

Lw - -lad + G)Wp'> = TT^^ = /<(« + gPp)»' + /^l ^ — jw + /u(a + gPp)w. 

We also have from the boundary conditions satisfied by w that w{po) = and 

3 



WpiO) 



-A-''^gp{Q)w{Q) + A-^'^gp{Q)w{0) 



Multiplying (13.341 ) by vv and integrating yields 

(w, Lw) - id(w, (a + gpp)w). 
Likewise, multiplying ( 13.341 ) by w and integrating: 

3/r „ ^. 1 \ , , . . I ii+G\ 



{Lw, w) - -Alail + G)wp\ , w\ = ip.{a + gpp)w, w) + \A-^ — W. wj + (//(a + gpp)w. 



w). 



We may integrate by parts to evaluate the second term on the left-hand side in the equation 
directly above, 

-tIWI +G)wp\ ,w\ = -Ml +G)Wp,Wp\- -fl(l +G)wp 



On the other hand. 



(vv, Lw) - {Lw, vv) = I [ 
J Pa ^ 



■w{a^Wp]p + w{a^Wpip 



)dp 



Wnd VV - w„a VV 



Adding the last three equations gives: 

-(a(l + G)wp, Wp\ - -a(l + G)wp 



w 



= (/((fl + gpp)w, w) 
■l+G 



I yUl — jw, VV j + I VVp VV - WpVV jfl 



Now as G(0) = for all A, G{Q) - 0. Thus the boundary conditions evaluated at /? = are the 
following: 

a^ (vvpvv - vvpvv) + —WpW = /I ' I - -/I" ' gpw + /I" ' gpvvvv j 



-A^^\A-^^hp^^) + \A'l\A-^^h^^) 



= 0. 
Thus, dropping the boundary terms, we have 
3 



-\a{l + G)wp, Wp\ = {Ji{a + gpp)w, w) + \A-^ — W, w\ 

But, a (and a + gpp) are strictly positive for admissible A. Hence, when ju < 0, we must have 
ju > as claimed. This completes the lemma. D 
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Applying these lemmas, we have, ultimately, the following result: 

Lemma 3.9. (Location of A*) Under the hypotheses of the previous lemmas, the solution 
A* of/j-iA") — —I (whose existence is given by Lemma [3741 ) is unique. Moreover, A* < Aq. 

Proof. By continuity of fx and the preceding lemmas, there exists some A* with fiiA*) - 
-1. Moreover, since /i is monotonically increasing when ^(A) is in any sufficiently small 
neighborhood of -1, this A* must be unique. 

We now prove that the size conditions given are sufficient to ensure that Aq + /!,. If 
A() < -2Bniin + eo, then we are done. Likewise, if eo = 0, then we are in the constant density 
case and this result is already known. So we shall assume /Iq > eo > 0. 

Fixing A - Ao,let<p e S' he given by 

f 1 + G(r; Aq) , 
Jpt, (Ao + G(r;Ao)y'^ 

Then, using (|3.27t we estimate 



Now, 



(loT^a-^-^- 



m'-' ' ''' -'"^ -''-'"'' 



But, for /l(), we know from (|3.16t that Y(po) = j^n- Combining these estimates we see that 
the numerator of the Rayleigh quotient 'R{-;A()) is dominated by 



•Jpn 



-gp(0)4>(0f + \ Yp(pr'<Pp(pfdp < -^1^ + 2Y(po) = 0. 
ho 2Y(po) 



Thus ix(Aq) < 0. 

Now let (f) e S minimize 'R(<p; /lo). Then multiplying the equation satisfied by (p and 
integrating by parts we find that, for po < p <0, 

rp rp 

a^cpcpp = -n(Ao) (a + gpp)4>^dr + a^c^ldr. 

Jpo Jpo 

Therefore, since we have already shown that -i-i{Ao) < 0, we see that a^(p(pp is a positive and 
increasing function of p. 

If we instead multiply the equation satisfied by (phy (p(\ + G)"' and integrate, we arrive 
at the following identity 



-gp{0)c^{Of+ r -^cf>ldp- f ^^cf>4>pdp = niAo) f ^ 

J Pa i + Cj Jn„ (I +Lrr Jn„ I + 



—d>^d 

G)^ Jpo l+G^ ^- 



Note that we have used the fact that G(0) = 0. To show that yi{A^)) > -1, therefore, we need 
to prove that 

-gpiQ)c^iQf+ c^^dp- —^M,4p+ —^4,^dp>0. (3.35) 

Jp„ I + (j Jpo (I + uy Jpo I + Lr 
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First we observe that, for any e S^ , 

\2 



= ( I (pp(p)dp\ 



1 r fl^^ 2 



Thus the first two terms in (13.351 1 give a nonnegative contribution. 

Note, however, that the third term is nonpositive, since Gp > 0. As we have seen, a^4"Pp 
is increasing, therefore we can bound the third term from below by 

Jpo (1 +Gr Jpo (1 +G)^ 

> -4gp(0)<p(Of\po\\\Gp\L (3.36) 

> -^gWppUpoWOf. (3.37) 

On the other hand, integrating the fourth term on the left-hand side of ( 13.35b by parts gives 
—^<P^dp = A{0)<p{Of - 2 —-rHpdp + --^^'dp, 

Jp„ i + G Jpg i + G JpQ (i + G) 

where Aip) :- J ia(r, Aq) + gpp{r))dr. Note that the final term on the right-hand side above 
is nonnegative. 

The next task is to estimate these quantities. In that regard we note 

-2 ^^dp>-A\pQ\A{0)ef\ (3.38) 

Jp„ 1 -I- G 

Then, again exploiting the fact that a^(p(pp is increasing, we use ( I3.37b -( l3.38b to find 

- r " ^' 4>^pdp+ f ^^^:^4>^dp>{\po\-'A(0)-4\\pp\\-4A(0)e^"^gp(0))\po\,p(Of. 

•J Pa (1 + G) Jp„ 1 -H G 

So long as the quantity in parenthesis on the right-hand side above is nonnegative, we have 
jL((/lo) > -1. But, since 

A{Q)>{e]^^-g\\Pp\U)\pa\>'^-\po\el'\ 

we have that 

Ipor'A(O) - WppU - 4A(0)eo'^V(0) > \\Pir'A{0) - \\pp\U > 0, 

where the last two inequalities follow from the definition of eo in ( 11.14b . D 
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3.4. Proof of Local Bifurcation. All that is left for us now is to verify the hypotheses 
of the Crandall-Rabinowitz bifurcation theorem presented in Q- As in the previous sections 
let the transformed fluid domain be 

R :- {{q,p) ■.Q<q<2n,pQ<p<Q}, 

with boundaries 

T := {{q,p) eR-.p^Q}, B := {{q,p) e R : p ^ p^}, 

and define 

X := {h € Cl:,"(R) -.h^OonB}, F = F; x 72 := C^.^W x €^^,"(7). 

Let h{q, p) :- H{p) + w(q, p). Then by the full height equation, w must satisfy the following 
PDE 

(1 + ^lq)(Hpp + Wpp) + Wqq{Hp + Wpf - 2WqWpq(Hp + W p) 

-g(H + w- d(H) - diw))(Hp + wp)'pp + (Hp + WpM-p) = in /?, (3.39) 

l+wl + (Hp + Wpf{2gp{H + w) - = on r, (3.40) 

together with periodicity in q and vanishing on B. 

Conforming to the framework of |7|, we introduce a nonlinear operator 

r = {TuTi) : (-2B™in + e^,^)xX^Y 

defined, for w eX, A> -2B„iin + eo, by 

T\{A, W) := (1 + w]q){Hpp + Wpp) + Wqq(Hp + Wpf - 2WqWpq(Hp + Wp) 

-g(H + w- d(H) - d(w))(Hp + wpfpp + {Hp + wpfpi-p) (3.41) 
TiiA, w):^\+w] + {Hp + Wp)\2gp{H + w) - Q). (3.42) 

Note that by definition of the laminar solution H{-;A), we have !F(/i, 0) = 0. For future 
reference, we evaluate the Frechet derivatives ;Fim , Tiw at w = 0: 

Ti^ ^dl+ Hid] + 3Hll5{-p)dp - 3g(H - d(H))Hlppdp - gHlpp{\ - d) (3.43) 
T2„ = {2gpHl + 2Hp(2gpH - Q)dp^ 

(3.44) 



= (2gpA-'-2A'/^dp^ 



In order to show that the eigenvalue at A* is simple we must characterize the null space 
and range of the linearized operator 'Fh,(A*,Q). This is accomplished in the following two 
lemmas. 

Lemma 3.10. (Null Space) The null space o/;7\v(/l*, 0) is one- dimensional. 
Proof. Fix A - A*. We have, by our work in Lemma [34l that M(p) cos q e N(Tw(A*, 0)), 
where 7V() denotes the null space. We need only prove uniqueness to complete the lemma. 



26 SAMUEL WALSH 

Let m in the null space be given. Then, since m is even and 27r-periodic in q, we may 
expand 



m(q, P) = 2^ cos (kq)mk(p). 



k=0 

The fact that "Fwi^*, 0)m = then implies 



oo 

^ {dl + Hid] + 3Hll3(-p)dp - 3g(H - d(H))Hlppd„ - gHlpp(\ - J)) cos (kq)m,(p) = 



k=(l 

in R, and 



CO 

^ i^gpHl + 2Hp(2gpH - 05^) cos (kq)mt(p)^: 



= 0. 

T 



k=Q 

We conclude that each term in the series above must vanish. But, note that 

d(mk(p) cos (kq)) - mk(0) T cos (kq)dq - 0, 
Jo 

hence, for each A: > 0, 

{-gpp + Hp^k^ + dp{Hfdp)]mk = 0, pQ<p <Q. 

From the second series and the boundary conditions satisfied by m, we have that for A: > 0: 

dpnik - gpA^^^^nik -0, on p -0, 

along with 

mil - 0, on p - Pa. 

For k - \, then, m must be a constant multiple of M, as both satisfy the identical ODE 
boundary value problem. Moreover, if A > 2, we see that nik = 0. Otherwise, the Rayleigh 
quotient in the minimization problem will be -k^ < -1, which contradicts the choice of A* . 

The difficulty arises in the case k - 0, since the (/-term persists. Using the fact dimo) - 
mo(0), we have that mo satisfies 

{-gPp + dp(Hp^dp)}mQ = -gppnioiO), pQ< p <Q 

along with the same boundary conditions as before on /? = and p - po- From this we wish 
to conclude that mo = 0. To see why this must be the case we rewrite it in the form: 



(fl^m[))' = gp'iniQ - mo(0)) = gp m^,' 

Jo 



{r)dr 

where p e {pQ, 0) and prime denotes differentiation with respect to p. Integrating this equation 
and using the boundary condition at /? = we find 

a^niQ = gp(0)mo(0) + I gp'(r) I mQ(s)dsdr 
Jo Jo 

= gp(0)mo(0)+ I ( I gp'(r)drjm'Q(s)ds 

C" 3 / 

= gp(0)mo(0) + I k{p, s)a(s) m(^{s)ds 

Jo 
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where the kernel, k, is given by 

k(p, s) :- I gp'(r)a(s) dr, po < p < s <0. 

This is nothing but a Voherra-integral equation of the form 

cf>(t) = C+ I kit, s)(f>(s)ds, pQ<t<s<0 (3.45) 

Jo 

with (p{Q) - C. As everything here is smooth and we are working on a compact interval, 
standard theory of integral equations gives a unique solution for - (p(p; C). Fix a value of 
C. Then multiplying (13.45b by a^^ and integrating we find 

npo nl'o nl'o r-t 

I <^(f; C)a{ty^dt = C I a{ty^dt + I a{ty^ I k{t, s)^(s; Qdsdt 
Jo Jo Jo Jo 

rpo rpo 

= C I a(ty^dt + I K(t)(f)(t; C)dt, 
Jo Jo 



where 



■-[ 



K{t) :- I k(s,t)a(s) ^ds. 



Now we note that we must have that a^m'^^ = 0(-; mo(0)), where (p is as above. Then, by 
the linearity of 0(-; C) in C, we see that 

fl^m;)(p) = gp(0)mQ(0)4>(p; 1), po < ;? < 0. 

But observe that a(pY^ - (A + G{p)Y^^^ and, recalling our argument of the previous section, 
we have that (p{p; 1) = 1 + G{p), as the kernels of the integral equations ( 13.281 1 and ( |3.45t 
(scaled so that C - 1) are identical. Hence integrating the above equation we find 

mo(0) = gp(0)mo(0) f -i^^^^fp = 27(0). 

Jp„ {A + G{p)yi^ 

Assuming oto(O) ^ 0, we arrive at a contradiction, since by the above identity 

Q{A) = 1 - 2gp{0)m) = 0. 



We conclude that m()(0) ^ only when A - Aq. However, Lemma I3T9I assures us that /I, < Aq, 
so this is impossible. 

We have, therefore, shown that mo(0) must be 0. This is an immediate consequence of 
the fact that m(0) = ot'(O) = while m satisfies a second-order, linear ODE. 

In summary, we have showed that all the modes k + \ vanish, thus miq, p) - m\{p) cos q. 
By Lemma 13741 see that ;«] is a constant multiple of M. Thus the null space is one-diminsional 
with generator M{p) cos q. D 



Lemma 3.11. (Range) The pair (^,!B) belongs to the range of the linear operator 
;7%v(/l*, 0) ijfit satisfies the orthogonality condition: 



I I J[a^<p*dqdp + - I 



Sa^0* = 
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where 0* generates the nulls space !?\,,(/l*, 0). 

Proof. Suppose first that (^, S) e F is in the range of T^^.{A*, 0). Then there exists some 
veX such that J{ = T\JA*, 0)v, and S = 5^2»v(/l*, 0)v. It follows that 

I I y[fl (f>*dqdq - \\ {{a Vp)p + aVqq - a gppv)(p*dqdp 

= I I {ia^fp)p + afqq - a^8Pp<f>*)vdqdp + I (fl\;,0* - a'^vfp)dq 

(a VpCp* - a v(p*p)dq, 



x< 



as on B, both v and 0* vanish. Now, on T we have a^ = /l^^^. Moreover, the ODE satisfied by 
(p* gives gpcp* = ^^^4'*n on 7". Thus, 

^ ^ {A^^^fpV-A'l\pf)dq. 

Combining this expression with the last gives necessity of the orthogonality condition. 

To demonstrate sufficiency we let any pair {J[, S) € F be given satisfying the orthogo- 
nality condition. We must show that there exists a solution m e X to the following: 

(-gPp(u - d(u)) + a^^Uqq + a^^{a^Up) - Ji mR, 
2(gpa-^u - aup) = S on T, (3.46) 

M = on B, 

Our methodology here will be to freeze the operator d, replacing it with a fixed real number. 
The resulting problem will be an elliptic PDE, so that we may use Schauder theory to extract 
solutions. 

We shall do this successively in several stages. First consider the problem, 

■ -ev(^> + -gppv^'^ + aM^q + fl-3(fl\'f ),, = y[ mR, 
2{gpa-\^'^ - avf) = S on T, (3.47) 

v*"^) =0 on B 

where v*"^^ is periodic in q. By standard elliptic theory, for each e > 0, problem (13.471 1 will 
have a unique solution. 

We claim, moreover, that these solutions are bounded in C^l"{R). By contradiction sup- 
pose otherwise. Then there is some sequence of e — > such that l|v^^^||(;i+a(m -^ oo. For each 
e, put M*'^^ := v*'^Vlk'^'llci+"(R)- Thus m''^* has unit norm for all e. Given that v'*^' solves (I3.47l i. 
we have additionally 



and 



-gppu'^^ + a-\fq + a-\a\f)^ -^ in C^t^iR) 



gpa-\<'> - auf ^ inC;er(r). 



Observe that ( 13.471 ) is uniformly elliptic. Applying Schauder estimates ensures that the se- 

-'2+Q' 

^per 



quence {m^'^*) is uniformly bounded in Cp^"{R). So by compactness we have a subsequence 
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converging strongly to some u e C^^^{R). By an argument identical to the previous lemma, 
we can show that u is in the null space of "FwiA*, 0), and hence a multiple of cp*. This follows 
from the fact that problem immediately above is T'uv without the stratification term. As we 
have seen, when we expand any element in the null space, this term will drop out in all modes 
A: > 1 . Similarly, for the k = Q term, we may simply bypass the integral equation argument 
(as we are, so-to-speak, already given that the right-hand side of the interior equation is zero), 
and use the same Rayleigh quotient manipulations to conclude uo = 0. So indeed, m is a 
constant multiple of <p*. 

But by definition of v*^*^' we have, 

ff Jkr'fdqdp = ^^ a^fi - ev*"' + TUA*, OW)dqdp. 

Applying the exact same manipulations on the second term in the integrand above as we did 
in proving necessity yields 

= — e II a v^'(p*dqdp — \\ a u ^' (p* dqdp ^=> = I I a u(p*dqdp. 
JJr JJr JJr 

However, we have found that u is in the linear span of 0*, so the equation above implies u = 0. 
This contradicts the fact that ||m||^i+„,^) = 1. Hence {v*'^^) is bounded. 

It now follows that {v^^^) has a strongly convergent subsequence in CL.^{R). The Umit, 
denote it v, will satisfy ( 13.47b with e = in the sense of distributions. But, again appealing to 
standard elliptic regularity theory, this implies that v e X. 

We have thus shown that for each (y[, S) satisfying the orthogonality condition, we may 
find a unique solution in X to the following problem: 

{-gPpU + a^'^Uqq + a^^{a^Up) - J{. mR, 
2{gpa-^u - aup) = S on T, (3.48) 

M = on B. 

Now to return to equation ( 13.461 ) we make the following observation: if the pair {JM, S) 
satisfies the orthogonality condition, then so does (^ - f(p),!B), where / is any smooth 
function. To see why this is the case recall that we have proved that 0* is of the form (p*iq, p) = 
M(p) cos q. Therefore 

I I (^- f)a^4>*dqdp = I I J[a^fdqdp - \ f{p)a{pf M{p)dp \ cos qdq 
JJr JJr J Pa Jo 






^a^ ^*dqdp 



Sa-" 



so, indeed, the orthogonality condition holds for (^ - f{p), S). 

Fix cr e R. In light of our previous observation, and our proof of the uniqueness and 
existence of solutions to ( I3.48l l in X, it follows that there exists a unique solution u^°'^ of 



( 



-gPpU^-> + a-^u%''+a-\a'uf\ 


= yi- 


- gPpO- 


inR, 


2(gpfl-2„(^) _ au^p) = S 






on T 


u^-^^ = 






onB 



(3.49) 



We may therefore define a mapping ^F : R — > R by ^(cr) := diu^"'^), where «*""* solves ( |3.49t . 
Suppose that ^' has a fixed point cr. Then we may simply substitute d{u^°'^) for cr in ( |3.49l l 
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to find that m^°"' is the sought after solution to (|3.46l l. Thus it suflices to prove that ^' has a 
(unique) fixed point. 

Let cr, T be given with cr + t and let m*°"\ m*'"* solve ( 13.491 1 for cr, t respectively. Then 
subtracting the equations satisfied by m*°^* and m*'^' and dividing by cr - r we arrive a solution 
to the following problem: 

-gPpV + a^h'qq + a'^ (fl^Vp) ^ = -gpp in R, 

2 [gpa-^v - avp) = on T, (3-50) 

V = on B. 

By applying the Schauder estimates for Dirchlet and oblique boundary conditions, moreover, 
we have that the solution, v, of this equation is in X. 

By contradiction, assume that d(v) - 1 . Then, by virtue of the fact v solves (13.501) . we 
have that v e NCFwiA*, 0)). It follows from our previous lemma that v is a constant multiple 
of 0*. But this is a contradiction, since 

d((p*) = M(0) f cos qdq = 0. 
Jo 

Thus we are assured d(v) + 1 . 

Fix cr, h eR. Examining (13.491 1, we readily see 

"Via- + h)- mcr) f / m^'^^'*' - m^'"' \ 

— h — ^ i \ — h — j^^ ^ ^^^^- 

Hence *P is differentiable everywhere and *I" = d(v). Therefore the function cr i-> ^{cr) - cr 
is monotonic with derivative uniformly equal to d{v) - 1 ^t 0. So for some cr* we have 
- i/'(cr*) - cr* . This proves the existence of a fixed point of ^, and thereby a solution to 
( 13.46b . The lemma follows. D 

Finally, we must ensure that the so-called transversahty or crossing condition holds. 



Lemma 3.12. (Technical Condition) If (p* generates the null space of 'F(A*,0), then 

Proof. First we calculate the mixed Frechet derivatives of'F?AA-A*,w- 0: 
TxA«{^\Q) = -(1 + G)a-^d] - 3(1 + G)a-^ pi- p)d p - 3gYa-%dp 
+ 3gY(l + G)a-%dp + lg(i+ G)a-'pp(l - d) 



TimM', 0) = ( - 2gpA-^ - A-'l^dp"^ 



By the previous lemma, it suffices to show that the pair {Tiji^iA* ,Q)(p* ,'F2Aw(A'' ,0)(p*) does 
not satisfy the orthogonality condition. Equivalently, if we put 



(1 + G)a-'{fY - 3(1 + G)a-'p{-p)4,\p 

3 



^gYappffp + 3gY{l + G)a-'ppf<f>; + -g(l + G)a-%(ffyqdp 
I ( - 2gpa'^{4>*f - -af4>*p\dq. 



+ 
It 
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our lemma will be proven if we can show H ^ (again, because d{(p*) - 0). We will 
demonstrate that, in fact, H < 0. To keep our notation concise, let H = Si + . . . + H7, where 
E, denotes the i-th term in the sum above. It is clear a priori that Hi > 0, while H5, He < 0. 

Consider the boundary terms He and H7. On T, we have a - A^^^, and (p* satisfies 
gp(f>* = A^^^cp* . Therefore, 

Sv = -]^^ A'l^ff^dq = A^A-'gpiffdq = ^Hfi. (3.51) 

Now, by the ODE satisfied by 0* in R, we have 

(a + gpp)f = {cr'4>*p)p = 3a(J3i-p) - gYpp)fp + a^(ppp, 
whence, 

3a'^ gYpp(/)*p = a^*j,p - (gpp + a)ar^(p' + 3a'^f5(-p)(f>*p. 
Substituting this expression into H4 yields the following, 

3gY(l + G)a^^pp(f>*^*pdqdp 

(a{\ + G)0V;^ - {gpp + a)a-\l + G){ffyqdp - H2 

J (a{\ + G)ffpp - gPpa-\\ + G){ffYqdp - Hi - H2. (3.52) 

Next, consider the quantity 

H2 + H3 + H4 = [T (fl(l + GW^*pp - gppa-\l + G){ff 

-3gYapp4>'(f>*\dqdp - Hi . (3.53) 

Calculating that Gp - -2gYpp, and Op = a^^(J5{-p) - gYpp), we integrate the first term by 
parts to find: 



-/I^ 



H2 + H3 + H4 = JJ {(a-'gY(l + G)pp - gaYpp)fcf,; - a-%-p)(l + G)4>*cf,; 

- fl(l + G)((f>*pf - ga-^(l + G)pp((f>*)Adqdp + i a(l + G)f<p*pdq - Hi 
J^ ( - ga-\l + G)pp(ff -a(l+ G)(<p;fyqdp 



+ -(H2 + H3+H4)-Hi-2H7, (3.54) 

where we have used the familiar fact that G(0) - 0. 

Now, simply combining identities (13.51b . (13.521 1. (13.53b . and (13.54b gives 

„„„„„„ 5„ 

.n = ill + £12 + ^^3 + -4 + ■^S + --6 

= Hi + ^( ^ ( - ga-\l + G)pp(rf -a(l+ GXf^fyqdp - Hi - 2H7) + H5 + ^Hg 
" ~2^' ~ 2 11 ""-^ "^ GXfpfdqdp + -He. 
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Thus H < 0, and the claim is proven. D 

We are now ready to prove the local bifurcation theorem. 



Proof of Theorem 15. 7] By definition of the laminar solutions, //(■; A), we have TiA., 0) = 
for all A > -2Btnin + eo- Moreover, by the regularity assumptions on p and/3 as well as the 
definition of X, we see that T^, Tk, 'Faw and ;Fuvv exist and are continuous. Applying Lemma 
[OOl and Lemma [37m moreover, we see that N(F„{A*, 0)) and Y \ 'R{'F,M\ 0) are both one- 
dimensional with the former generated by 0*. Finally, Lemma |3T2l shows T^^iA* , 0)0* is not 
contained in ■??(5^v('i*,0)). 

Thus we have satisfied all the hypotheses of the Crandall-Rabinowitz theorem on bifur- 
cation from a simple eigenvalue (cf. fT\ Theorem 2). This allows us to conclude that there 
exists a C' local bifurcation curve 

C'Q^{A{s),w{s):\s\<e), 

for e > sufliciently small, such that (_A{Q), w(0)) = {A*, 0) and 

{{A, w) e 1/ : w * 0, T{A, w) = 0) = Co 

where 1/ is some neighborhood of {A* , 0) e [-2Binin + eo, °°) x ^- Moreover, for \s\ < e, 

w(s) = s(f>* + o(s), inX. (3.55) 

Furthermore, since h - H + w and Hp > in R, we may restrict our attention to a smaller 
C' -curve C[g^ c C'q containing (A*, 0) and along which hp > in R. Then by Lemma ITTl we 
have the existence of a C'-curve, Cioc, of solutions to (ll.lb - (ll.6l l n 

4. Global Bifurcation. We now prove that we can continue the curve Cj'^^ whose ex- 
istence was established in Theorem 13.11 Continuing with our notation from the previous 
section, denote 

X:^{he Clt"(R) : /i = on B), F = Fi x Y2 - C^^TW x Clt"(T), 

where the subscript "per" indicates 27r-periodicity and evenness in q, R is the rectangle 
(0, In) X (po, 0) and T = (0, In) x {/? = 0). Next define 

^ = (^i,^2):RxX-^y 

by 

gm := (1 + hpipp + hg^hl - Ih^hphp^ - g(h - d(h))h]fip + hlpi-p) (4.1) 



g2{Q,h):^(\+h] + hl{2gph-Q)^ 



(4.2) 

p=0 



Recall that here t/ : X — > R denotes the linear operator mapping an element of X to its average 
value on T. We remark, also, that the laminar flow solutions {H(A),X) found in Lemma [3721 
satisfy Q(H{A\ Q{A)) = 0, for all A > -2B„in + eo- 

In order for our arguments to have any traction we will need to make strong use of a 
priori estimates. It will often prove convenient to consider uniformly elliptic differential 
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operators that approximate Q, or more specifically, Qi . For any cr e R we therefore define 

g^'^^ : R X X ^ F by 

g'^pih) := (1 + h])hpp + h,,hl - 2h,hphp, - g(h - cr)hlpp + hl/3(-p) (4.3) 



gf\Q, h) := GiiQ, h)^(\+hl + hlilgph - 0); 



(4.4) 

p=0 



That is, we replace the li-term of Q\ with the real number cr. As will seen later, in the function 
space we shall work in, this defines a uniformly elliptic diff'erential operator for each cr e R. 
This is essential, for it allows us to exploit Schauder theory in order to prove the compactness 
properties we need. 

Let 5 > be given. In order to ensure the uniformity of the ellipticity of ^j and 
obliqueness of @^ we will work in the set 

Os := 1(2, h)e^xX:hp> 5mR, h< ^-^ on r|. 

Likewise, we put 

Sg-.^ closure in R x X of {(Q, h)eOs: g(Q, h) = 0, h^ * 0), 

and let C^ be the component of Ss that contains the point {Q*,H*), where Q* := Q{A*), 
H* :- //(■; A*). Thus C'^ contains the local curve C[^^. 

For later reference we compute the Frechet derivatives of ^ and ff^'^K 

§ih{h) = ihqhppd, + (1 + hl)dl + 2hph,,dp + hid] 

-2(hpghqdp + hpqhpdq + hphqdpdq) 

-3gpp(h - d{h))hldp - gpphlil -d) + 3hlp(-p)dp (4.5) 



&2h{QJi) = (2/1,5, + 2gphl + 2hp(2gph - Q)d,)j. 



(4.6) 



GTh^h) = e^hih) + 3g(d(h) - O-Wpppdp + gpph'pd 
= 2hghppdq + (1 + h])dl + 2hphgqdp + hid] 
-2(hpghgdp + hpqhpdq + hphqdpdq) 

-3gpp{h - cr)hldp - gpphl + 3h]fi{-p)dp (4.7) 

0th^Q^h)^g2k(Q,h). (4.8) 

In the spirit of Rabinowitz 1.20.1 . we shall apply a degree theoretic argument to get a global 
continuation theorem in the form of an alternative result: 

Theorem 4.1. (Global Bifurcation) Let 6 > be given. One of the following must hold: 
(i) C'g is unbounded in R x X. 

(ii) C'g contains another trivial point {Q{A), H(A)) € 7~, with A + X'. 
(iii) C^ contains a point (Q, h) € dOs- 

Observe, however, that Q is not a compact perturbation of identity, which rules out using 
the classical Leray-Schauder degree (see, e.g. |,11J). In light of the nonlinear boundary 
operator 02, we instead employ a variant degree theory developed by Healey and Simpson 
(cf. lfT4l ). In order to do so, we must first establish two lemmas on the topological properties 
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of the map Q. The structure of the arguments in both cases will be to begin by using elliptic 
estimates on ^^°"' (which is uniformly elliptic in Og). Then, taking advantage of the fact that 
the operator d can be easily estimated in X, we reinsert the stratification term and are able to 
make conclusions about Q. 

We emphasize that these lemmas are key. Once Q has been shown to be admissible in 
the sense of Healey-Simpson degree (cf. Definition 4. 10 of fT45), the proof of Theorem 4. 1 
is identical to that of the homogeneous case in |5|. Indeed, the technical obstacles presented 
by the stratification term are completely confined to the following three proofs. 

Lemma 4.2. (Proper Map) Suppose K is a compact subset of Y and D is a closed, bounded 
set in Os, then @^^(K) C\ D is compact in R X X 

Proof. First we prove that for each Q, cr e R, the operator h i-> 0'^'^\Q, h) is uniformly 
elliptic and oblique in Og. The former follows from the fact that the coefficients of the higher 
order terms in Q\h satisfy: 

A{\+h])hl-Ah]hl^Ahl>6\ 



Notice that the bound here does not in any way depend on cr or Q. Similarly, we have that 

^2 (/i) is uniformly oblique as the coefficient ofhp 



the boundary operator h i-> g^ W is uniformly oblique as the coefficient of hp satisfies 



\{2gph - Q)hp\ >6^ on r 

for(e,/i)€^. 

Let {ifj,gj)} be a convergent sequence in F = Fi x F2 and assume that for each j > 1, 
(fj,gj) = @{Qj,hj) for some (Qj,hj) e O's with [hj) bounded in C^^^iR), Qj bounded in R. 
We wish to show that there exists a subsequence of {{Qj, hj)] convergent in R x X. 

Towards that end we denote 6j := dqhj, for j > I- Then, differentiating the relation 
between {Qj, hj) and (fj, gj), we find by ( 14.11 ) 

dqfj ^dgSiihj) 

= (1 + {d,hjf)d],ej + (5,(1 + {dghj)^))dlhj + {dphjfdlOj + {dlhj)d,{dphjf 
-2(d,hj)(dphj)dpdg6j - 2(dlhj)(dphj)(dpd,hj) - 2{dqhj)(dpd„hjf 
-gidphjfppOj - g(hj - d{hj))d,{dphjf + d,{dphjfii{-p), in R, V j > 1 . 

Grouping terms above we may rewrite this in the form 

(1 + (d,hjf)dlej + {dphjfdl6j - 2(d,hj)(dphj)dgdp0j = 

djj + F(d,hj, dpdqhj, d^hj, dphj, d^hj, d(hj),pp), in R (4.9) 

where F is the cubic polynomial dictated by the previous expression. Notice that we have 
picked up a dependence on d(hj) and pp that was not there in the constant density case pre- 
sented in 151. This will not be an issue, however, as these terms can still be bounded in 
X. 

We may likewise differentiate the equation on T to discover that, for each j > 1, 

dcigj^dqg2(Qj,hj) 

= 2(dghj)d,9j + 2(dphj)(2gphj - Qj)dpej + {dphjf(2gpd,hj), on T. 

Equivalently, 

2idghj)d,9j + 2idphj)(2gphj - Qj)dpej = d,gj + G(d,hj, dphj,p), on T, V7 > 1 (4.10) 



STRATIFIED WATER WAVES 35 

where G is the quadratic polynomial determined by the previous equation. Finally we observe 
that 

Oj^Q, onB. (4.11) 

Since {ifj,gj)} is convergent in Cp^"{R) x Cp^"{T) we have that the sequences {dqfj}, 
{dqgj} are Cauchy in C"(R) and C}^^"(T) respectively. Moreover, by the boundedness of {hj} 
in Cp^"(R) (and thereby the boundedness of d{hj) in this space) we have that F is uniformly 
bounded in C}^^"{R), and G is uniformly bounded in Cp^"(T). It follows that each of these, 
viewed as a sequence in j, is Cauchy. Then, by the compactness of the embeddings, the right- 
hand sides of equations ( 14. 9t and (14.10b are pre-compact in C"(R) and C},^"{T), respectively. 
Possibly passing to a subsequence, we may take both to be convergent in these spaces. 

We now consider differences 6j - 6k, for j,k > 1. By ( 14.9b we have 

( 1 + {d,hjf)dl(ej -ek) + (dphjfdliOj - Ok) - 2(d,hj){dphj)d,dp{ej - e,) = Fj^, in r 

where by our arguments in the previous paragraph we know Fjit — > in C"(/?). Similarly, 
from ( 14.11b have that 6j - 6k vanishes on the bottom and on the top (14.10b tells us 

2{dqhj)d,{ej - Ou) + 2(dphj)(2gphj - Qj)dp{9j - Ok) = Gju, on T. 

Here G jk — > in C'^"(r), again by the considerations of the preceding paragraph. We now 
apply the mixed-boundary condition Schauder estimates to the differences 9j - 6k to deduce 
that \\6j - 0k\\c-+arj\ — » as J, fe — > oo. We have shown, therefore, that all third derivatives of 
hj are Cauchy, except possibly for dphj. To demonstrate the same holds for {dp2j}, we use 
the PDE to express d^hj in terms of the derivatives of order less than or equal to two: 

dlhj = (1 + {d,hjfr'i^fj - d]hj{dphjf + 2(dghj)(dphj)(dpdghj) 

+g(hj - (T){dphjfpp + (dphjfpi-p)^, inR, j>l. 

But we have seen that the right-hand side is Cauchy in C'^"(/?), hence {d^phj] is also Cauchy 
in C"(R). We conclude that the original sequence, {hj}, had a convergent subsequence in 
Cp^"(/J), and hence {{Qj, hj)} had a convergent subsequence in R x X D 

Lemma 4.3. (Fredholm Map) For each {Q, h) e Os, the linearized operator QhiQ, h) is a 
Fredholm map of index from X to Y. 

Proof. In the previous lemma we established that g'*""' was uniformly elliptic and oblique 
V(2, h) G Os and Vcr e R. We remark also that these bounds are independent of cr. 

Now let 1^ e Cj^+"(R) be given and fixcr,Q eR, h e X. Put 

0*" := dg {&t\Q, h)W) - (5,^[f (2, h)) W, for / = 1, 2. 

Here, by (dg^r^ (Q, h))[ij/] we mean differentiating the coefficients of ^[^ in q, then applying 
the resulting operator to if/. Then dytf' satisfies: 

ff%\h)d,i// ^ cp'-'^ onR 



g%\Q,h)dgif, = ^^^^ onT 
dqif/ = on B 
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which is a uniformly elhptic PDE with an obhque boundary condition. The classical Schauder 
estimates ensure the existence of a constant C > 0, independent of xf/, such that 

< mW(R) + \\d,{g'^^{h)n\c-(R) + \\d,{@''2t!{Q,h)n\o..(Ty 

On the other hand, we may express d^pif/ via the partial differential equation to arrive at an 
estimate for dtij/ of the same type. Combining we get 

C||(Allc3«(i) < ll«Allc2-(«) + ll^9^^"llc"(«) + Il5?<^*''llc-(r)- 

If we now apply the classical Schauder estimates to (/», we find that for some C > 0, 
independent of if/. 



C||«Allc2«(«) ^ ll-Allc-w + II^4M1*W'A)|Ic«(«) + \\dc,{Q^2h{Q,h)ilj)\\c. 



'(T)- 



Together with the previous line, this implies that there exists a constant C - C(cr) > such 
that, for all t// e C^^"(R) even and periodic in q that vanish on B, 

c(cT)\mx < ii-Aiic-w + iiMa'W'Aiif, + \\&^^^Q,h)nY,. (4.12) 

Of course, this is not quite the estimate we are after, since there is a lingering dependence on 
cr. In order to remedy this we make the following elementary, but very useful, observation: 
for any t// as above we have 



Wm < j Wq < ||iAllc«(«)- 



We are therefore able to estimate terms involving d easily in spaces of Holder continuous 
functions. In particular, 

< C||«Allc....(«)' (4.13) 

where the constant C above depends only on pp, \\h\\x and our choice of cr. Likewise, an 
identical argument gives the estimates 

\\di&ih{h)t// - 5;^i'^'(/i)iAllc°(R) ^ C\W\c-+«(Ry i^ P, q 

where C depends again on pp, \\h\\x and the choice of cr. Of course we do not need to make 
such arguments to estimate the §21, (^) term, as it identical to @2h{h). 

Combining these observations with our first estimate (14.121 ). we find that, for some C > 
and all 1^ e X, 

cm\x < ii<Aiic-"(«) + w&ihwnY, + \\g2h{Q,h)nY2- (4.14) 

The key point here is that X - Cp^"{R) is compactly embedded in the spaces who appear 
on the right-hand side of ( 14.141 ). To show that the null space of @h(Q, h) is finite dimensional, 
for example, it suffices to show that the unit ball is compact. But for any sequence {i/r,,) c 
{i/f G N{0h{Q,h)) '■ \\^\\ci+«(R) ^ 1). we have, by the compactness of the embedding, there 
exists a convergent subsequence [ijJnt] in C^^^iR). This subsequence is Cauchy in C^^"(R), so 
applying ( 14.141 ). we see it is also Cauchy in X. Finally, by completeness, {i/r,,, ) is convergent 
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in X. It follows that the unit ball in the null space is compact, hence the null space is finite 
dimensional. 

The proof that the range is closed is, likewise, standard. Again, the argument hinges on 
the compact embedding of C^^^iR) cc C^^^iR). For brevity we omit the details. 

Now, the Fredholm index has a discrete range, hence by the connectedness of Os it must 
be constant on this set. But QhiQ* , H*), which we denoted 'Fm.(/1*, 0) in the previous section, 
was shown to have a one-dimensional null space in Lemma 13.101 Moreover, by the orthog- 
onality condition of Lemma 13.111 its range has codimension one. Since {Q* ,H*) e Os by 
construction, the Fredholm index is uniformly along the continuum C'g. D 

Finally we prove a result characterizing the spectrum of the operator Q. 

Lemma 4.4. (Spectral Properties) (i) V5 > 0, 3ci , C2 > such that for all {Q,h) e Os with 
\Q\ + \\h\\x < M, for all ip &X and for all real ji > C2, we have 

cilltAllz </"^ll(A-//)(Ally, +fi'^\\Btlx\\Y, 

where A = A{Q, h) = ^ihih) and B = B{Q, h) = QuiQ, h). 

(ii) Define the spectrum S = 2(2, h) by: 

l.(Q, h) :— {/J. e C : A — ij. is not isomorphic from {i// e X : Btji - onT] onto Y\]. 

Then 2 consists entirely of the eigenvalues of finite multiplicity with no finite accumulation 
points. Furthermore, there is a neighborhood N of [0, +oo) in the complex plane such that 
2(/l, w) C\ N is a finite set. 

(iii) For all (Q, h) e Os, the boundary operator @2h{Q, h)from X ^> Y2 is onto. 

Proof. In this proof we follow the well-tread path laid forth by Agmon in |[T1. The 
only novelty here is in (i), where we must do a little work before we can apply the elliptic 
estimates. The remaining parts are completely standard; their proofs rely only on (i). Lemma 
I4.2l and Lemma 143] (cf. lU, ||5l or lfT4l . for example). For that reason we omit them here and 
devote our attention to (i). 

Fix cr e R, and let (A € X and// € C be given. Put 9 :- argyU, and suppose for some e > 0, 
l^l < 7r/2 + e. Consider the operator D^"^' := A^"''^ + e'^d^ onRxR, where A^'"' := gf^(h). 
Then, D*°"' is elliptic for each cr, and the boundary condition on the top is oblique, hence the 
complementing condition holds. Let f ; R — > R be a cutoff" function supported compactly in 
the interval / := (-1, 1). Put 

e(t) := e'l'^l"''^(0, 4>(t,q,p) := e(t)ilf(q,p). 

We apply the Schauder estimates in R^ with the boundary operator B to (p(t, q, p) to deduce 

C||0|lc3«(/xR) < II<^IIC"(/X«) + IIO^'Vllc-x/xR) + l|B0llc^«(/xr) (4.15) 

for some constant C > independent of i/r. 

A quick calculation readily confirms the existence of C, C > 0, depending only on (, 
with 

c^l^ < Ikllcnn) < c'pT, c^i^ < lkllc.«(R) < CV^. (4.16) 
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Using this, we can unpack ( 14.151 ) to derive the following set of estimates: 

iP*"Viic'«(/x«) = ||e(f)^*';>w[«A] -Mo<A + (i/^r^'r + r)e'vw"'v|| 



IC'"(/xfi) 



for lyul sufficiently large. We can estimate the second term in parenthesis by appealing to 
( 14.13b in the previous lemma, yielding 

IP^'Vllc'-x/xs) < CIjuI"? (||(^ia(/2) -;«)<AIIc'™(«) + ll<Alb«w), (4.17) 

for \fj.\ sufficiently large. 

Similarly, analyzing the boundary terms we find for some C > 0, independent of i^, and 
|yu| sufficiently large, 

\m\y^«(ixT) = Mt)g2hiQ,h)my«(ixT) 

<C\^i\'^my..^ry (4.18) 

On the other hand, using ( 14.161 ) and the fact C^^" cc C*^", for A; < 3, we can show that for 
some C > 0, independent of t//, and for all \jj\ sufficiently large, 

my.^m^Cfiimy..(Ky (4. 19) 

Now, inserting ( 14.171 ), ( |4.18t and ( 14.191 1 into ( |4.15l l, we arrive at the inequality in (i). D 

Let 'W be an open bounded subset of X and fix Q e R. By virtue of Lemma |42}Lemma 
14.41 we conclude that 0{Q, ■) : 'W — > F is admissible in the sense of [ 14 1. This enables us to 
use the generalization of the Leray-Schauder degree introduced in that paper, which we now 
briefly recapitulate. 

Let y e F \ g(Q, d'W) be a regular value of g(Q, ■). We define the (Healey-Simpson) 
degree of 0(Q, ■) at y with respect to "W by 

deg(^(e,.),nV,3'):= ^ (-l)''^"'', 
weg-'(iy}) 

where v(w) is the number of positive real eigenvalues counted by multiplicity of ^u, (2, w), 
subject to the boundary condition §2w(Q, w) = on T. By Lemma |4~41 v(w) is finite. Likewise, 
the properness of 0{Q, ■) established in Lemma |4~2l implies ^"'({y)) n 'W is a finite set. The 
degree is therefore well-defined at proper values. In the usual way, this definition is extended 
to critical values via the Sard-Smale theorem (this is valid by Lemma l43T l. 

Our interest in degree theory stems from the fact that the degree is invariant under any 
homotopy that respects the boundaries of 'W in the following sense. Let 1/ c [0, 1] x 'W be 
open. Define 1/, := {w e OV : (t, w) e ^^} and likewise ^^l, ■.^{we'W : it, w) e 5'ZY). 

Lemma 4.5. (Homotopy Invariance) The degree is invariant under admissible homo- 
topies. That is, suppose 1/ c [0, 1] X 'W is open. IfH € C^(1J; Y) is proper and, for each 
t € [0, 1], 'Hit, .) is admissible (i.e., satisfies the conclusions of Lemma I4.21 Lemma I4.4I ), then 
we say fi is an admissible homotopy and have 

deg (7^(0, ■),'i/o,3') = deg CK(1, ■),'l/i,3;) 

provided y i Wt, dV,)for all t e [0, 1]. 
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From our previous analysis it follows that 0{Q\, ■) i--> QiQi, ■) is an admissible homotopy. 
Thus, in light of Lemma |431 the degree will remain constant as we move along the contin- 
uum. Assuming that all three alternatives of Theorem |4T| fail, we use this feature to generate 
a contradiction. At this stage, we have reestablished all the relevant properties of ^ that were 
true in the constant density case. Repeating (verbatim) the proof in f5\ we obtain Theorem 



5. Nodal Pattern. We now seek to investigate the second possibility of the global bi- 
furcation theorem, namely that the continuum arcs back and intersects T at some point aside 
from A* . This will be done by assuming there exists some other laminar solution on C'^ and 
analyzing the nodal pattern to show it to be, in fact, equal to A* . In particular, we will be 
concerned with the vanishing of hq. We therefore work in the set Q := (0,7r) x {po,0) c R. 
Denote 



do., 
dQh 
dQ.r 
dill 



= {(q,Po) :g€(0,:T)}, 
= {i7T,p): peipo,Q)}, 
^{(0,p):pe(pQ,0)}- 



It follows that h - on dQ.b for all iQ,h) e C^, and for h e X periodicity and evenness 
in q yield hq = on dO.,- U dQ.i. In the following lemmas we will attempt to prove that for 
heCg\{{Q\H*)}v/e have additionally 



hq<0 

hqp < 
hqq <0 
hqq >0 



in Q U do., 
on 5Qi 
on dQ.1 
on dQ.r 



(5.1) 



and on the bottom corners of Q: 



hqqp{0,Po)<Q, hqqpin, pu) > 0, 



at the top corners, we have 



hqq(7:,Q)>0, V(0,0)<0, 



(5.2) 



(5.3) 



These inequalities define an open set in X. 
Our first result is the following: 



Lemma 5.L Properties (I5.1b -( l53] l hold in a small neighborhood of {Q*,H*) in R X 
Cp^"(Q) along the bifurcation curve CJ^^. \ {{Q* ,H*')] originating from the point {Q* , H*). 

The proof of this lemma does not make any special reference to the form of the operator 
Q, only the basic properties of the eigenfunction w* and the local bifurcation curve. Not 
surprisingly, therefore, it follow with virtually no modification from the homogeneous case 
treated in ||5]. We therefore omit it and concentrate on the remaining lemmas, which will 
require more finesse. 
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Differentiating the relation §\ (li) - by ^, for any choice of cr e R, yields fflv(h)[hq] = 
for (Q,h) e C'^, where we recall that 

0u(.h)[(P] = 2hghpp4,g + (1 + hl)<f>pp + 2hphggCf>p + /Z^,^,, 

-2(hpc,hq(pp + hqhp(pq + hphqCppq) 

-^gPpih - o-)hl(pp - gpphl(p + 3hl/3{-p)<pp. 

Likewise, differentiating the boundary relation we find that QihiQ, h)[hg] = 0, for (Q, h) e C^, 
where 

@2h{Q, h)[<l>] := 2hg(Pg + 2gphl<p + 2hp(2gph - Q),Pp. 



An important observation is that ^.T — though elliptic for given cr e R and h e (9, 



s ■ 



has a zero-th order term. We have dictated that pp < 0, so that the sign will go precisely 
"the wrong way", in the sense that the maximum principle will not hold in general. The 
argument will be saved by the observation that the (p under consideration is known a priori to 
be non-positive. If we denote 

'H^^\h):^gf^(h) + gpphl, 



then 'H^°'Hh) is uniformly elliptic with no zero-th order terms. Moreover, every non-positive 



1 that is a solution relative to Q^"^{h) is a subsolution relative to 'H^"'\h). 



Lemma 5.2. Properties ( I5.1l )-( l5.3l l hold along C'j \ {(Q*,H*)} unless 3A + X* with 
(Q(A),H(A))eC',. 

Proof. By the local bifurcation theorem we know that in a sufficiently small neighbor- 
hood of iQ*,H*) the bifurcation curve C'^ consists entirely of the curve C'^^. Suppose by 
contradiction that the statement of the lemma is false. As the previous lemma implies that 
(EB-dSJl will hold near (Q\ H*) in R x C^^+"(Q), we have the existence of some (Q, h) e C^ 
with hq ^ where at least one of the properties fail, although they hold on a sequence 
(2y> hj) — > {Q, h). Continuity then implies that Z;^ < on Q, as this holds for each hj. Since 
C^ c Os, we have that 0ff^\Q, h) is uniformly elliptic where we take cr :— d(h). 

We wish to apply exploit Hopf 's lemma to conclude the last three inequalities of ( 15. Il l 
hold. To do so, we first note that, as h^ < in O and h^ = on dQ.i U dO., U 5Q/,, we have 
sup^hq = 0. But then, by our comments above, we see that cp :- hq is a subsolution relative 
to the linear elhptic operator 'K'^^^/z). Applying Hopf's lemma to 'H^°'H}i), we conclude 



50 

5^ 



>0, 

(90 .Po) 



where (^o, Po) is any point where 4> vanishes, and v is the outward unit normal at {qo, po) (cf, 
for example. Theorem 2.15 of \\\2X ). In particular, since dy - -dq on cJQ/, and hq = on this 
set, we have hqq < on 3Q/. An identical argument shows hqq > on dQ.r, and hqp < on 

To establish the strict inequality hq < on 50,, we appeal to the nonlinear boundary 
condition, (I4.2l i. We already have that /z^ < in this region by continuity. Suppose that for 
some q() e (0, n) we have hq(qo, 0) = 0. Letting - hq we use our expression for^2/i(2' h)[(p]. 
Furthermore, 

2hp(2gph - Q)hqp ^ at(qo,0). 
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Again, by Hopf s lemma, we have that hqpiq^, 0) > 0. Moreover, as {Q, h) e O^, we are 
guaranteed hp > S > Q. As p > 0, this implies 2gph - Q onT . But then, the nonlinear 
boundary condition for the original problem 0{Q, h) = would imply I + h^ - 0, which is a 
contradiction. Hence we have the strict inequality and therefore ( 15. Il l holds for h. 

In order to produce a contradiction, therefore, all that remains is to verify the corner prop- 
erties ( I5.2l i and ( I5.3l l. By continuity we have that hqqp(n,po) > 0. Suppose that it vanishes. 
We have that h(q, po) - for all q, hence hq{n, po) - hqq{n, po) = hqqq{n, pa) - ... =0. 
Also, by evenness and periodicity in q we have hq{n, p) - for all p. Hence hqp(n, p) - 
hqppin, p) - . . . - Q, for all p. But then if v = Vqq + Vpp, where q - (\,Q), p - (0, 1), is any 
vector exiting Q. at {n, po), then 



dv 
and 



Vqhqq + Vphpp = 0, at (tt, po) 



= at(7T,po) 



5v2 

by the previous analysis. This violates Serrin's edge lemma producing a contradiction (see, 
for example. Theorem E.8 of 1121 ). An identical argument applied to (0, po) proves that ( 15. 2t 
holds. 

For ( 15.31 1 we consider the point (0,0). By evenness we know that hq{0, p) - Q for all 
p e (po;0). Thus hqp{Q,p) = for all such p. By continuity, moreover, we are guaranteed 
that hqq{0, 0) < 0. By contradiction suppose that hqq{0, 0) = 0. Then differentiating the 
relation QihiQ, h)[h\ - twice yields 



hqhqq + gphqhp + hphpq(2gph - g) = 0, On dn, 



and 



f^lg + '^g'hg + gphgghl + Igphqhpq + hlg(2gph - Q) 

+hphpqq(2gph - Q) + 2gphqhphpq = 0, on dD.,. 

Evaluating both of these at (0, 0) and taking our assumption into account, we have all terms 
involving hq, hqq and hpq drop out. Hence 

hphpqq{2gph-Q)^0, at (0,0). 

But we already have seen that 2gph - Q cannot be zero, as this leads to a violation of the of 
the nonlinear boundary condition. Moreover, hp is bounded uniformly away from in Os, 
so we must conclude that hpqq{Q, 0) = 0. As before this violates Serrin's edge lemma, as any 
vector V leaving Q. through (0, 0) can be written as v = Vpp + Vqq, and therefore 



dv 
and 



dhq 

^- = Vqhqq + Vphpq = 0, at (0, 0) 



5v2 



= at (0,0), 



since hpqq{Q, 0) = 0. This proves that ( 15.3b holds for h, hence we have a contradiction estab- 
lishing the lemma. D 
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Combining the two previous lemmas we see that the nodal properties (|5.1| l-( l53] l will 
hold along the continuum C'g unless it intersects the laminar curve somewhere other than A*. 
We now conclude this section by showing this cannot occur for A to the right of -2Bniin + £o ■ 

Lemma 5.3. If a trivial solution (Q(A), H(A)) € C'g, with A > — 2fi,„,„ + eo, then A — A*. 

Proof. Let (Q{A), H{A)) e C'g be given. Clearly if A - A* we are done, so suppose 
that this is not the case. As we know locally to (Q*,H*) that C'g = C[^, we may assume 
without loss that there exists a curve C'g(A, A*) c C^ originating at (Q*, H*) and terminating at 
{Q, H) that does not intersect 7~ except at the points A, A*. Let {{Q{A„), h„)] be a sequence in 
C'f,{A, A*) \ T such that {An, h„) -^ (A, H(A)) in R x C^+"(Q.). Thus dgh„ i 0, and g(Q, h„) = 0, 
Vn > 1 . We claim that dqhn < for n > 1 . By the first lemma of this section this holds 
true for all h in some sufficiently small neighborhood of {Q* , H*) in X. The arguments of 
the previous lemma, however, can be directly applied to show that the nodal properties hold 
along C'^iA, A*) \ {(Q{A), H{A))). That is, if we assume that for some h e C^(^, A*), (l5Jli-(l53Tl 
fail, then we may again choose a sequence of functions on C'l^iA, A*) approaching h where 
the properties hold. The previous lemma can then be applied without further modification to 
produce a contradiction. This proves the claim. Thus dqhn < 0, for n> \. 

We now differentiate the relation 0{Qn, h„)[hn] = with respect to q to find that, for each 
n > 1, 

= dqffi (h„)h„ 

= (1 + idqh„f)dldqh„ + 2{dqhn){d]h„){dlhn) + {d\h„){d ph„f 

+2(dlh„)(dph„)(dpdqh„) - 2(dlh„)(dph„)(dpdqh„) - 2(dqh„)(dpdqh„f 
-2(dqhn)(dphn)(dpd^ghn) - 3g(hn - d{hn))(dphnf{dqdphn)Pp - g(dqhn){dph„fpp 

+3idph„fidqdphnm-p) 

where v„ :- \\dqh„\\^l^^ -dqh„ and o-„ := d(h„). Likewise, we may diff^erentiate the boundary 

relation to conclude that QX^ (Qn, h,^\yr\ - 02h{Qn^ h„)[vn\ — 0, for all « > L 

For brevity denote (g, H) :- {Q(A), H{A)) and cr := d{H). We observe that since h„ -> H 
in €^^"(0.), we must have cr,, — > cr. Combining this with the convergence of /?„ — » // in the 
C3+"-norm, we see {gff;-\h„)} is a Cauchy sequence in £(€^^"(0), €"(0.)). Then for /, ; > 1, 
we have 

m%'\hd - ^i?\^))v;llc.(n) < K,\hd - M?'(^)llx(c-"(n).c"(n)) ^ 0, as /, j ^ ex,, 

L Likewise, {g'^'^"'(2„, h„)} is a Cauchy sequence in £(€^^"{1), C'+"(r)), 
so that 

mf^\Qi, hi) - g'^'\Qj, hj))vj\y..^T^ ^ 0, as /, ;■ ^ ex,. 

By linearity we have 

@T'\Qi, hdlvi - vj] + (g[^^\Qi, hi) - G'^'\Qj, hj))[vj] = 0, V/, g>\. 

Notice that for each /, j > 1, we have that Qp\Qi, hi) is uniformly elliptic, hence we apply 
the Schauder estimates to the above equation and conclude that there exists a generic constant 
C > independent of /, j satisfying: 

C\\vi - V;||c2+„(n) < ||v; - V;||c(.(n) + \\{@t'\Qi,hi) -Q^^' (2>/ij))v;ilc"(n)xc'+"(r)- 
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Our previous estimates imply that the last term on the right-hand side vanishes as /, y — > oo. 
On the other hand, as the {vj} are uniformly bounded in C^^"(Q.), the compact embedding 
of C^^"{Q.) cc C*'(Q) guarantees there is subsequence convergent in C°(Q). But then the 
inequality above implies this subsequence converges in €^^"(0.) as well. By abuse of notation 
we identity the subsequence with {vj} itself. We note also that as d(vj) - d{dqhj) - 0, 
for all j > I, that the limit must also have mean zero. Then we may let m be given with 
dqm e C^^"(Q) and Vj — > dqm in this space. It follows that ||5^m||c2+„(j^) - 1, and 

Qt\Q^H)[dqm] = Qh{Q,H)[dqm\ = 0, 

where the first equality follows from the fact that dqin has mean zero so the d term in Qh 
vanishes. 

As we have argued before, evenness and periodicity in q together imply that dqhj vanishes 
on dD.r U dQ.1 for j >\. Moreover, as hj = on dQ.b, we have that vy = on dQ.b U dQ.i U dQ.r. 
We have, additionally, that dqhj < in Q, for each j > 1 . Together these imply that dqin = 
on 5Q/, U 5Q/ U dQ.r and dqin < in Q. But dqin has unit norm and satisfies the uniformly 
elliptic equation Qj'^ [dqm] - 0, hence we may apply the maximum principle to conclude 
dqin < in O. 

Mirroring our derivation of the eigenfunction at (Q*, H"), we expand dqin in a sine series 

oo 

dqiniq, p) = 2_^ fu(p) sin kq, q,p e R. 

k=0 

As 0hiQ,H)[dqin] - 0, each mode in the series expansion must also satisfy this relation. 
Evaluating QhiQ, H) we may drop all terms involving q derivatives to find 

§ih{Q,H) ^dl+ Hjdl - 3gpp(H - d(H))Hldp - gppHl + 3Hlp(-p)dp, 

02h(Q, H) = {igpHl + 2Hp(2gpH - Q)dp) = l^gA' - A^'^dp^ . 

In particular, for k - I find 

= (^2 + Hjdl - 3gpp(H - d(H))Hldp - gppHl + 3Hlf5(-p)dp)fi , Vp e (po, 0), 
or in self-adjoint form: 

dp(Hfdpfi) = (H;,' + gpp)fi , ^p e (po, 0). 
Similarly, applying 02hiQ, H) we find 

/i(/7o) = 0, /;(0) = gp{o)A-'i\fm- 

Returning to our original notation, this implies 

-l^nSuA)>p{A) 

where we recall H is the Rayleigh quotient defined in (13.231 1. If p(A) < -1, then -1 is not 
the minimum eigenvalue, so f\ is an excited state of the above Sturm-Liouville problem. We 
conclude by standard Sturm-Liouville theory that f\ vanishes at some point in (pQ, 0). But we 
know from explicit computation that 



fiip)^- (dqin)(q, p) sin qdq<0, pe{po,0) 

J^ Jo 

as dqin < 0. This contradicts the existence of nodes, so we must have instead that p(A) - 
-1 - piA"). By Lemma [T91 it follows that A - A* us desired. D 
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6. Uniform Regularity. In the previous section we considered the second alternative of 
Theorem 14.11 Now we wish to delve further into the first. We make the central aim of this 
section the establishment of uniform bounds in the Holder norm along the continuum C'g. The 
ultimate product of our efforts will be the following theorem. 

Theorem 6.1. (Uniform Regularity) Let 6 > be given. If sup (j^ q.^^^^, snpj hp and 
suP(/,,0eq 121 are finite, then supy, gj^q II'^IIc^+^cr) is finite. 

To treat the second derivatives we shall employ a suite of a priori estimates for nonlinear 
elliptic equations with oblique boundary conditions due to Gilbarg and Trudinger (cf. ifTTl ). 
We shall not, however, need these results in their full generality, so instead we streamline 
them into a single statement. 

First note that due to the periodicity in q, the domain R can be considered as embedded 
on a torus, allowing us to ignore the seeming lack of smoothness at the corner points q - 
0, 2ji. Next, in keeping with Lieberman and Trudinger's notational framework we consider a 
differential operator F = F{z, ^, r) e C^(R x R^ x S, R) and boundary operator G = G(z, f) e 
C^(R X R^, R). Here S denotes the space of 2 x 2 real symmetric matrices, D denotes the 
gradient operator and D^ the Hessian. Both of these are taken as acting on the space of 
smooth real-valued functions on R which are 27r-periodic in the first variable. As usual, we 
say that F is elhptic at (h,^, r) G R x R^ x S provided that the matrix Fr := [dF/drij]i<ij<2 
is positive definite at that point. Moreover, if Ai and A2 denote the minimum and maximum 
eigenvalue of Fr, respectively, then F is said to be uniformly elliptic provided that the ratio 
A2/A1 is bounded. The boundary operator G is said to be oblique at a point (q, 0) if the 
normal derivative ;t' := Gf ■ (0, -1) is positive at (^,0) for all (h,^) e R x R^. 

Assume that F and G are given as above and consider the corresponding nonlinear elliptic 
boundary value problem 



(6.1) 



Theorem 6.2. (Lieberman-Trudinger) Let h e C^(R) be a solution to (16. Il l that is In- 
periodic in the first variable and for some constant K > satisfies \h\ + \Dh\ < K in R. Suppose 
further that for some positive constant M the functions F — F(z, ^, r) and G — G(z, ^) satisfy 
the following structural conditions: 

A2<MAi, (6.2) 

|F|, IFfI, \F,\ < MAi(\r\ + 1), (6.3) 

Frr < 0, (6.4) 

|G|, |G,|, IGfl, |G,,|, |G,f|, IGffI < Mx, (6.5) 

for all (z, ^, r) € R X R^ X § such that \z\ + \^\ < K. Then there are positive constants p — p{M) 
and C — C{K, M) such that h e CperiR) and 

This is a nontransparent restatement of the original result presented in I.17J , so we take a 
brief aside to outUne its development. The general idea will be to incrementally estimate the 
Holder norms, beginning at CP^^^ and working upwards to C^^^. For clarity, all citations from 



F(h, Dh, D^h) -- 


= 0, inR 


G(h, Dh) = 


on/7 = 


/i = 


onp = po 
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IfTTll are marked with the abbreviation L-T. 

Proof. Let F e C^CR x R^ x S, R) and G e C'^iWx R^R) be given and consider the 
boundary value problem ( 16. 11 1. First we wish to show that assuming ( I6.2l )-( l6.5b . F and G 
meet the hypotheses of L-T Theorem 2.1, namely structural requirements (L-T Fl), (L-T 
F2), and (L-T G2). But observe that conditions dO and dOJ imply (L-T Fl) and (L-T F2) 
respectively, where we taking fiQ to be a constant. Moreover, ( 16.5b ensures that condition 
(L-T G2) holds, again taking yUo to be a constant function. So indeed, the hypotheses of L-T 
Theorem 2.1 hold. Applying the theorem, we have for any C^iR) n C^iR) solution h with 
\h\ < M, there exists positive constants // - jj(n, M,fj.()) and C = C(n, M,ij.q,R) with (16.21) 



[h^R < C. 

Here [■]ft;R denotes the Holder seminorm. If we slightly abuse notation and drop the depen- 
dencies on n and R, it follows immediately that there exists a constant C = C(K, M) such 
that 

Next we wish to show that structural requirements (L-T 4. 1 )-(L-T 4.4) are a proper subset 
of (IE2b-(l63Tl. Clearly, dO i is equivalent to (L-T 4.1). Also, taking 6-0 and noting that 
F^ = 0, we see that (IgU l gives both (L-T 4.2) and (L-T 4.3). Finally, dOI ) along with the 
obliqueness of G implies (L-T 4.4), with an appropriate choice of /ii . This verifies all of the 
hypotheses of L-T Theorem 4. 1 ; hence for some //' > and C > 0, both depending on M 
and K, we have 

[Dh]^.,R < C. 

By abuse of notation we can therefore, choose ^ and C positive constants such that 

Again we note that these constants depend only on K and M. 

For the higher derivatives we shall need to make use of L-T Theorem 5.4 and L-T Theo- 
rem 6.2. Conditions (L-T 5.1) and (L-T 5.2) coincide exactly with conditions (L-T 4.1) and 
(L-T 4.2), respectively. Similarly, condition (L-T 5.4) follows from the positive definiteness 
of F and (L-T 5.22) from ( 16.5b . We have already seen that (L-T 4.3) holds, so applying the 
second statement of L-T Theorem 5.4, there exists a positive constant — which we again 
denote C — depending only on K and M with 

sup \D^h\ < C. 

R 

Lastly, in order to apply L-T Theorem 6.2 and get an estimate of the full C^^'^-norm we need 
only verify that sti-uctural conditions (L-T 6.1), (L-T 6.2), (L-T 6.3) and (L-T 6.4)' hold. But 
(L-T 5.1) and (L-T 4.3) together imply (L-T 6.1), and (L-T 6.4)' clearly follows from (L-T 
5.22). Positive definiteness of F, moreover, gives (L-T 6.2) for a sufficiently large fi2. As 
this implies we have a bound on the C^-norm of A^^F, (L-T 6.3) also holds for large enough 
//3. Applying L-T Theorem 6.2, then, and combining with the previous estimates yields the 
theorem. D 

That settled, we return to the proof of the uniform regularity theorem. 
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Proof of Theorem \6.1\ We begin by proving the L'^{R) norm of h is uniformly bounded 
along the continuum. To do so we note that by (|2.9l l 

< h{q, p)=y + d< 7]iO) + d, 

and thus the definition of po in ( ll.lOl i yields 

J](0) + d< ° = \po\ sup hp 

mf^ ^Jp{c - u)) ^ 

where the last equality follows from (12.10b . Thus assuming the hypothesis of the theorem, 
that is that hp is uniformly bounded in L°° along the continuum, then the above inequalities 
together give the boundedness of the L'^-norm of h. 

Next, to bound hq we fix Q and differentiate in q the full height equation (|2.17t to find 

- 3hpghl(S(-p) = (1 + h^g)hppq + Ihqhqqhpp + Ihqqhpqhp + hqq^h], 
~^'^qq'lp'^pq ~ ^'^q'^nq ~ ^'^q'^p'^pqq 

-ghqhlpp - 3g(h - d)hpqh],pp, in R, (6.6) 

Ihghqq + 2hpghp(2gph -Q) + Igph^h] = 0, on /7 = 0, (6.7) 

hq-Q, onp - p(). (6.8) 

Consider now the function k{q,p) :- h(q,p) - d + eqe"P, where e > and « e N will be 
specified later. Inserting k - eqe'^P in place of h - d in ( I6.6I 1 we have 

O(e^) = -6eqne"Pkpgkp/3(-p) + 3kpgkl/3(-p) - 3eqne"Pkl/3i-p) + Ik^k^^kpp - 2en^qe""kgkg, 

— €6 KppKqq + (i + tC^)lCppq — Z€€ K-qKppq — (i + tC^)6n C — ZlCqlCp^ + Z€e Kpq 
+H€n€ rCpqlCq — ZlCptCqtCpqq + Z€C i<-pi<-pqq + Z^Uqe rCqlCpqq + tCplCqqq — ZEUqe f^pf^qqq 

-gkqlc'pPp + 3geqne"''kqklpp + gee"''k],pp - 3gkkpqk^pPp + 6geqne"''kkpkpqPp 
+3g€qe"%qklpp + 3gene"Pkklpp, inR. (6.9) 

At an interior maximum for kq in R we would necessarily have 

f^pq ~ '^qq ~ ^5 ^PP — ^' Q1^ — ^' f^qppf^qqq — ^qqp- (O.iU) 

Thus evaluating ( 16. 9t at such a maximum point gives the following simplified equation: 

O(e') = kppq + kppqikq - ee""f - Ikpqqikq - ee"P){kp - enqe"P) 
+kqqq(kp - enqe"Pf - gkqlkjpp - 3eqne"Pklpp] 
-ee"P [r? + r?k] + 3/3(-p)nkl - gklpp - 3gnkklpp]. (6.11) 



Likewise, inserting k into ( I6.8I 1 we find kq - ee"''" on p = pQ. In particular this implies that at 
the interior maximum kq is positive. We may therefore choose n € N to ensure that the quanti- 
ties in square brackets above are positive. In fact, it suffices merely to have rp-+3l3{-p)nk^ > 
and 3£qne"'' - kp > in R to achieve this, taking into account the signs of pp, k, kp and kq. 
But combining this fact with the signs implied by ( 16.10b yields a contradiction, since for e 
small we must have the right-hand side of (16.11b is negative. We conclude kq does not have 
an interior maximum. On the top we have from ( 16.7b 

1 + (kq - ef + {kp - eqnf{2gp(k + d - eq) - Q) ^ 0, p^O. 
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By assumption h^ - (kp - enqe"PY is uniformly bounded along the continuum, and by our 
previous argument h-d - k- eqe"P is as well. But then, as the maximum of kg occurs on the 
boundary, and we have shown that on the boundary kg is controlled by these quantities, we 
must have that the maximum of kq is uniformly bounded along C'^. Clearly, the same must 
hold for /i^. We can repeat the argument but instead taking the auxiliary function A; = h-eqe"'', 
which will show that the minimum of hq is likewise bounded. 

For the second derivatives we fix cr e R and apply Theorem 16. 2 1 to Q^'^K To do this we set 

Fiz,^, r) := (1 + ^2)^22 + rnf2 ' ^^i^irn - giz - (r)f^Pp + ^ll5(-p) 
Giz, a := 1 + ^^ + f2(^gpz -Q) onp^Q, 

where ^ = (^1,^2) £ R^ and r - irn,ri2,r22) £ R^- By using cutoff functions. Theorem 
I6.2l is applicable in a subset in which solutions exists a priori and for each choice of cr. In 
particular, ^2 ^ <5 implies ( 16. 2t . ( 16. 3t is immediate and as Frr - 0, ( 16. 4t holds. Finally, as the 
only modification to G that occurs when considering the variable density case is the addition 
of the constant p(0), a trivially modification of the constant density argument in |5 1 shows that 
(16.5b holds. It follows that for each cr e R, THEOREM l6.2l applies to the solution of0^°'^[h] - 0. 
Of course, this introduces the possibility that the constant C will depend on cr. However, 
every h € C'^ is such a solution for cr = d(h), and recall that we have cr - d{h) < sup^ \h\. 
Thus we may select K independently of cr. Moreover, since in the definition of F, cr does 
not occur in any of the r^ terms, it does not affect Ai and A2. Indeed, the characteristic 
polynomial satisfied by the eigenvalues is 

A^-(l+f,+ ^2)A,. - (1 + ^2)^2 _ 4^2^2 ^Q . ^ J 2 

which does not depend on cr. This is an obvious consequence of the fact that the effects 
of density variation in the height equation are consigned to the lower order terms. Finally, 
the term containing cr in F can be easily estimated in terms of K. This fact, along with 
the uniform boundedness of the L°°-norms of h and Dh along the continuum, allows us to 
conclude that there exists C > such that 

sup Pllc^..(s) < C. 

Finally we establish third derivative bounds. Denoting 6 - hq, we can differentiate the 
height equation in q with fixed Q to find 



(6.12) 



with Cpgl'(R) coefficients and right-hand side / a Cper(^) function. Since we have already 
proved h is bounded uniformly in Cp^^iR), it follows that the right-hand side of (16.121 1 is 
bounded uniformly in Cper(^) along the continuum. As we have seen (in the proof of Lemma 
14.21 for example), so long as hp is bounded uniformly away from along the continuum, the 
problem above will be uniformly elliptic with a uniformly oblique boundary condition. We 
may therefore apply Schauder theory to obtain uniform a priori estimates of 6 in Cp^^(R) over 
all of Q. 

It remains only to bound hppp in Cpf-^iR). By means of the height equation ( 12.13b . we 
may express hpp in terms of the lower order derivatives of h: 

hpp = -(1 + h^y {h^pPi-p) + hqqhl - Ihqhphpq - g(h - d(h))hlpp) 



ii + hl)epp-2hphgepq + hlegq 

= f(h,hp,hq,hpp,hpq,hqq) 
hq6q + gphle + {2g(TK{h] - Q)hp9p 


in/? 
= on /:> = 


61 = 


on p = po 
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But the right-hand side above is in C^g^iR) by the arguments of the previous paragraphs. 
Thus hppp in C'^^iR) as desired. 

To transition back to the original Holder exponent, a, we merely note that if h e Cp^^(R), 
then it is certainly in C^^"(R). The arguments we have used to derive the third derivative 
bounds in no way relied on the particular value of//, so running through them with a instead 
we find h g C^+«(^). d 

7. Main Result. With the regularity established in the previous section, we are now 
prepared to begin the task of proving the main theorem. Theorem 1 1.1 1 Theorem 14. II gave us 
three possibilities for the continuation of the local bifurcation curve. Then, Theorem 6.1 of 
the previous section allows us to conclude that if maxj dph and Q remain bounded along the 
continuum, then h remains bounded in X. We now unpack the remaining alternatives and 
consider their significance in the context of the original problem (ll.7l )-( fT78T l. 



Proof of Theorem \l.l\ Fix (5 > 0. If C^ is unbounded in R x X, then at least one of the 
following must occur 

• there exists a sequence (g„, /;„) € C'^^ with lim„^co Qn - °°\ 

• there exists a sequence (g„, h„) € C|^ with lim„^co maxjdphn - co; or, 

• C'g contains another trivial point (Q(A), H(A)) e T, -IB^nin < A < -2Bjnm + fo < ^*, 
where the second alternative follows from Theorem 16. H and the third from the lemmas of §5. 
If C'j contains a point of dOg, then either 

• there exists a (Q, h) € C'j with dph = 6 somewhere in R, or 

• there exists a {Q, h) e C'^ with h - ^— somewhere on T. 

Thus, Theorem |4T| tell us that at least one of the five alternative above must occur. Consider 
the first alternative. By the definition of po in ( ll.lOl i we can estimate 

inf ^/p(x^(c - u(x,y)) < —^^ < sup ^Jp(x,y)(c - u(x,y)) 



and thus appealing to ( 12.51 ) we deduce 
Q ^ p(0, ?7(0))(c-m(0, 77(0))f +2gp(0, ?7(0))(?7(0)+qO < sup p(c-m)^+ ^8P(0,m)po ^^^^ 



If for some 6 > the first alternative holds, we conclude from ( 17.1b that for the corre- 
sponding sequence {u„, v„,p„, ri„) of solutions to ( ll.7b -( ll.8b . either sup^j— p„(c - u„) — > oo, or 
inf^— p„(c - Un) -> 0. 

Similarly, if the second alternative holds, then simply by rearranging the change of vari- 
ables in ( 12.121 ) we have 

dph„ = 



^|p^{c - u„) ' 

hence raax-^dph„ — > oo iff inf^j— p„(c - u„) -^ 0. 

Next suppose that for some decreasing sequence (J„ — > the fourth alternative holds. 
Again appealing to ( 12.12b we see immediately that the corresponding sequence (u„, v„,pn, rjn) 
must satisfy sup;^— p„(c - m„) — > oo as n — > oo. 

Finally, if we have a sequence 5„ — » for which the fifth alternative holds, then there 
exists a sequence {{Qn, h„)} with (Q„, h„) e C'^ for each n > 1 and supj{2gph„ - Q„) — > 0. We 
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have shown previously that dqh vanishes on dQ./. Thus if we evaluate the boundary condition 
at the crest (0, 0), we find 

Qn - 2gp„{Q)hn(Q, 0) < e„ - 2gpni0)hniq, 0), V<7 € [0, Itt]. 



idph„(0,0)) 

Therefore dph„{Q,0) — > oo as n — > oo. In particular, we have max-^dph„ — > oo so that the 
fourth alternative implies the second. 

By definition the family of continua C^ indexed by 5 > is increasing as 6 decreases. 
We may therefore define C :- [Js>oC'g to be the maximal continuum. Then by the consid- 
erations of the §2, and in particular Lemma 2. 1, there exists a connected set C of solutions to 
(ll.lb - (ll.6b corresponding to C. The arguments of the preceding paragraphs, moreover, imply 
that either supT^— m„ — > c, infT^— m„ — > -oo along some sequence in C, or C contains more 
than one distinct laminar flow. This completes the proof. n 

Remark. In the homogeneous case it was proved in ||5| that if inf^^ u„ — > -oo along 
some sequence in the continuum, then necessarily sup^ u„ ^ c along some sequence. In 
other words, there are waves on C whose speed is arbitrarily close to the speed of the wave 
profile. 

We conjecture that the same holds true in the heterogeneous case. Suppose, on the con- 
trary, that for some sequence infg— m„ — » -oo, but m„ remains uniformly bounded away from 
c on C. Let (Qn, h„) be the corresponding sequence of solutions to the height equation. Then 
we can show that sup-^,— 5y,/i„, sup^,— /i — > 0. In particular, this implies ||77n||oo, and d{h„) —> 
in the limit. As the fluid domain is bounded between these to values, this entails that, mov- 
ing along the continuum, we have regions of the fluid moving arbitrarily fast to left, while 
simultaneously the fluid is vertically pinching off. This is a consequence of the fact that, if 
u„ -^ -oo, then in order to keep a constant pseudo-mass flux po, the domain must be vanish- 
ing. Giving the pathology of the scenario, we strongly believe it does not occur, though we 
have been unable to rule it out. n 

Acknowledgments. The author wishes to thank W. Strauss for his indispensable guid- 
ance and support during the writing of this paper. 

REFERENCES 

[1] S. Agmon, On the eigenfunctions and on the eigenvalues of general elliptic boundary value problems, Comm. 

Pure Appl. Math., 15 (1962), pp. 1 19-147. 
[2] C.J. Amick, Semilinear elliptic eigenvalue problems on an infinite strip with an application to stratified fluids, 

Annali Scu. Norm. Sup. Pisa CI. Sci. 4, 11 (1984), pp. 441^99. 
[3] C.J. Amick and R.E.L Turner, A global theory of internal solitary waves in two-fluid systems. Trans. Amer 

Math. Soc, 298 (1986), pp. 431^84. 
[4] J. Bona, D.K. Bose and R.E.L. Turner, Finite amplitude steady waves in stratifled fluid, J. Math. Pares et 

Appl., 62 (1983), pp. 389^39. 
[5] A. CoNSTANTEM AND W. Strauss, Exact Steady periodic water waves with vorticity, Commun. Pure Appl. Math., 

57(2004), pp. 481-527 
[6] A. CoNSTANTiN AND W. Strauss, Rotationul steady water waves near stagnation, Phil. Trans. R. Soc, 365 

(2007), pp. 2227-2239. 
[7] M. Crandall and P. Rabinowitz, Bifurcations from simple eigenvalues, J. Func. Anal., 8 (1971), pp. 321-340. 
[8] B. Cushman-Roisin, Introduction to Geophysical Fluid Dynamics, Prentice-Hall, Englewood Cliffs, NJ, 1994. 
[9] M.L. Dubreil-Jacotin, Sur la determination rigoureuse des ondes permanentes periodiques d'ampleur flnie, 

J. Math. Pures et Appl., 13 (1934), pp. 217-291. 
[10] M.L. Dubreil-Jacotin, Sur les theoremes d'existence relatifs aux ondes permanentes periodiques a deux di- 
mensions dans les liquides heterogenes, J. Math. Pures et Appl. 9, 16 (1937), pp. 43-67. 
[1 1] I. FoNSECA and W. Gangbo, Degree Theory in Analysis and Applications, Clarendon Press, Oxford, 1995. 



50 SAMUEL WALSH 

[12] L.E. Fraenkel, An Introduction to Maximum Principles and Symmetry in Elliptic Problems, Cambridge Uni- 
versity Press, Cambridge, 2000. 
[13] D. GiLBARG AND N.S. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer-Herlag, 

Berlin, 2001. 
[14] T. Healey and H. Simpson, Global continuation in nonlinear elasticity. Arch. Rational Mech. Anal., 143 

(1998), pp. 1-28. 
[15] K. KiRCHGASSNER, Wavc-solutions of reversible systems and applications, J. Diff. Eq., 45 (1982), pp. 1 13-127. 
[16] K. Lankers AND G. Friesecke, Fast, large-amplitude solitary waves in the 2d Euler equations for stratified 

fluids. Nonlinear Anal., Theory, Methods and Appl., 29 (1997), pp. 1061-1078. 
[17] G. LiEBERMAN AND N.S. Trudinger, Nonlinear oblique boundary value problems for nonlinear elliptic equa- 
tions. Trans. Am. Math. Soc, 295 (1986), pp. 509-546. 
[18] R.R. Long, Some aspects of the flow of stratified fluids. Part I: A theoretical investigation, Tellus, 5 (1953), 

pp. 42-57. 
[19] A.S. Peter and J.J. Stoker, Solitary waves in liquids having nonconstant density, Comm. Pure Appl. Math., 

13(1960), pp. 115-164. 
[20] PH. Rabinowitz, Some global results for nonlinear eigenvalue problems, J. Func. Anal., 7 (1971), pp. 487- 

513. 
[21] G. Stokes, On the theory of oscillatory waves. Trans. Cambridge Phil. Soc, 8 (1847), pp. 441^55. 
[22] A.M. Ter-Krikorov, Theorie exacte des ondes longues stationnaires dans un liquide heterogene, J. 

Mecanique, 2 (1963), pp. 351-376. 
[23] R.E.L. Turner, Internal waves in fluids with rapidly varying density, Ann. Scuola Norm. Sup. Pisa CI. Sci. 8, 

4(1981), pp. 513-573. 
[24] R.E.L. Turner, Traveling waves in natural systems in Variational and Topological Methods in the Study of 

Nonlinear Phenomena, Progress in Nonlinear Differential Equations and Their Applications 49., V. Benci 

eds, Birkhauser, Boston Basel Berlin, 2002, pp. 114-131. 
[25] M. Yanowitch, Gravity waves in a heterogeneous incompressible fluid, Comm. Pure Appl. Math., 15 (1962), 

pp. 45-61. 
[26] C.S. YiH, Dynamics of Nonhomogeneous Fluids, Macmillan Company, New York, 1965. 



